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steady  BK>tion«  (3)  Discussion  of  the  problem  of  viscous 
fingering,  (4)  A  statement  of  the  physics  of  Immiscible 
fluid-fluid  displacement,  (5)  N€t%iork  model  theory,  (6)  Oenerallsed 
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PHBFACE 


Shortly  after  the  seventh  progress  report  f«as  coaq^leted* 
the  principal  investigator  (ehile  retaining  supervisory 
control)  withdrew  his  active  participation  in  the  Contract 
work  because,  even  with  a  suppleaental  allowance  received  in 
July,  there  %«as  not  enovtgh  money  to  cover  more  than  the  sumswr 
salaries  of  the  permanent  project  workers  (research  associates 
and  assistants) .  In  due  course  it  was  subsequently  learned 
that  the  support  of  the  work  would  terminate  Avtgust  31,  and 
the  task  to  prepare  a  final  report  followed  as  a  consequence. 

With  the  return  of  the  principle  investigator  to  the 
Illinois  campus  in  septesiber,  it  was  considered  premature  to 
prepare  the  required  final  report.  For  example,  the  Masters 
thesis  ««ork  of  N.  Chaudhari  was  to  be  inritten  in  the  fall 
as  based  on  the  experimental  data  obtained  during  the  sunaier. 

H.  Fare,  who  had  a  half-time  appointment  from  February 
through  August,  was  still  in  the  process  of  crystallizing  his 
ideas  about  network  model  studies,  integration  of  the  Stokes- 
Wavier  equations,  amd  about  the  meaning  of  Brittin*  s  %fork. 

The  situation  with  w.  Rose  %fas  that  he  had  a  number  of  partly- 
finished  manuscripts  which  needed  additional  work  before  they 
could  be  considered  for  publication.  In  fact,  the  only  project 
%K>rlcer  in  a  position  to  write  a  suimnary  statement  was  R. 
Channapragada  (id>o  was  committed  in  advance  to  tenninate 
at  the  end  of  the  siamner,  and  who  was  committed  therefore  to 
leave  a  statement  of  his  progress  over  the  two'^ear  period 
of  his  contract  work) .  unfortunately.  Or.  Channapragada 
did  not  find  the  time  to  prepare  the  «ranted  definitive  sumsiary 
of  his  work. 
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forn  of  th«  final  progress  report  as  attached  re¬ 
flects  on  the  factors  nentioned  above.  A  cosqpletely  rewritten 
Fara  and  Rose  paper  is  Included,  together  with  a  statsaent 
by  Fara  about  his  ideas  on  prograsming  solutions  of  the 
Stokes-Navier  equations,  solutions  of  certain  categories 
of  network  aodel  studies,  and  other  questions.  An  Abstract 
of  the  thesis  of  N.  Chaudhari  is  included.  A  nuaiber  of 
interim  statessmts  are  offered  by  W.  Rose,  and  finally  there 
appears  the  brief  statements  of  Or.  Channapragada. 

In  due  course,  supplemental  statements  will  be  offered 
to  complete  the  material  as  offered  in  this  final  report. 

That  is,  we  take  the  position  that  our  continuing  'work  under 
new  sponsorship  had  its  origin  in  the  CNL  contrsct  work,  and 
that  an  obligation  rests  on  us  to  provide  complete  statements 
at  the  end  of  all  projects  started  under  the  CML  sponsorship. 
Specifically,  the  final  version  of  Chaudhari' s  thesis  will  bs 
offered,  together  with  revisions  of  certain  of  the  manuscripts 
of  W.  Rose  and  others  contained  herein. 
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snwwwf  or  proghiss 

Ovmx  th«  tWD-y«ar  period,  Septeaiber  1959  through 
Auigust  1961,  a  project  titled  "Analytic  Study  of  the  BUhavior 
of  Flulda  in  Poroua  Hedia"  was  aupported  by  the  Chemical 
Warfare  Laboratories  (Directorate  of  Research)  at  the 
university  of  Illinois.  Ihis  work  is  identified  by  Contract 
Nusber  DA«18-108-405-ClfIi-517,  and  has  led  to  the  preparation 
of  eight  quarterly  progress  reports  of  udilch  this  is  the  last 
auid  final  report. 

During  the  sxibject  period,  some  20  students,  post«docto> 
rates,  consultants  and  faculty  were  employed  for  varying 
periods  of  tisie,  cumulatively  resulting  in  the  support  of 
some  seven  nantyears  of  effort.  Three  masters  thesis  %rere 
completed  under  contract  sponsorship,  and  some  one  dosen 
manuscripts  were  prepared  and  are  being  prepared  and  approved 
for  publication. 

It  is  clear  from  the  foregoing  that  the  contract  funds 
%wre  used  largely  for  support  of  students,  of  which  six  were 
doctorate  candidates.  The  principal  research  associate, 

R.  Channapragada,  for  exasqple,  received  his  PhD  degree 
(in  mathematics)  in  June  1961.  Less  than  two  percent  of  the 
contract  funds  has  been  used  for  the  purchase  of  equiptment 
and  laboratory  supplies. 

The  contract  work,  while  it  has  led  to  some  new  scientific 
findings  (all  of  %«hich  are  either  now  pid>lished  or  soon  to  be 
published),  has  been  particularly  meaningful  in  forming  a 
foundation  for  future  activities  asad  research  studies  to  be 
undertaken  in  this  Laboratory  forthwith  under  National 
Science  Foundation  and  other  Department  of  Defense  support. 

Wo  list  the  major  scientific  achievements  of  the  contract 
work,  as  follows: 


(55 


1.  Thm  developntnt  of  a  battar  undaratanding  of  the 
Intareonnactlon  batwaan  tha  nlcro  and  nacroseoplc  aapaeta 
of  tha  unsteady-atataa  of  alxtura  flow  In  poxoua  aadla. 

(of.  Bosa  and  Channapragada,  ZZZ  and  ZV*;  Boaa  VZZZ) ; 

2.  An  Invaatlgation  of  tha  xola  of  alectrokinatlea  In 
fluid  flow.  (cf.  straat,  V;  straat  aitd  Sao«  ZV; 

Stawart  and  Straat,  ZV;  Straat  and  Stawart,  ZV; 

Straat,  ZZZ,  ate.). 

3.  An  axparlaantal  invaatlgation  of  the  velocity  and 
acceleration  dependency  of  moving  contact  anglea. 

(Chaudhari,  VZZZ;  Roaa,  VZZ;  Haina,  VZZ;  Chaudhari,  VZZ; 

Hoaa  and  Chaudhari,  VZ;  Roaa  and  Haina,  V;  Roaa, 

Ramagopal  and  Haina,  ZV;  Haina,  ZV,  ate.). 

4.  An  invaatlgation  of  tha  meaning  and  uaa  of  tha 
Brittin  aquation  to  daacriba  capilliury  flow,  (cf .  Fara 
and  Roaa,  VZZZ;  Para  and  Roaa,  VZZ;  Tung  and  Orawa,  VZ; 

Roaa  and  Fara,  VZ;  Fara,  Roaa  and  Tung,  V;  Tung,  ZV;  ate.). 

5.  A  conaldaration  of  tha  solution  of  network  modal 
problama.  (cf.  Fara,  VZZZ;  Fara,  VZZ;  Tung  and  Roaa,  V; 

Fara  and  Roaa,  V:  Roaa,  ZV;  ate.). 

As  for  the  continuing  studiaa  in  prograaa  in  thia  laboratory 
which  rapraaant  an  outgrowth  of  tha  CMXf 517  %iork,  wa  refer 
particularly  to  tha  atatsmant  of  Roaa  (VZZ)  titled)  *1ha  Baaic 
Problem  Raf ormulatad" .  Bara  it  is  made  clear  that  to  even  begin 
to  meat  tha  broad  dbjactivaa  contained  in  tha  original  pra-> 
aabla  to  tha  CMZi-517  project,  an  appeal  auat  be  made  to  tha 
fundamental  fluid  maehanica  and  aurfaca  chamiatry  of  porous 
madia  ayataaai.  Indeed,  it  was  in  tha  spirit  of  this  view 
to  which  tha  various  general  contributions  of  Channapragada 
ware  directed  (cf.  contributions  in  VZZZ,  VZ,  V  and  ZV  in 
particular) . 

*Rafara  to  tha  contribution  of  Rosa  and  Channapragada  as  it 
appeared  in  tta  third  and  fourth  quarterly  progress  reports,  ate. 
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lf«  conolud*  that  tha  pzoblanMi  of  hlghast  priority  to 
which  wo  should  now  direct  our  attention,  include: 

1.  mnnerical  means  to  acconqplish  a  "brute-force” 
integration  of  the  equations  of  Bavier-Stokes,  with 
major  nonlinearities  ratatned  (cf .  Fara  and  Oeiger,  VIZZ) . 

2.  A  further  study  of  how  charged  double  layers  (i.e. 
the  electrokinetic  effects)  control  the  shapes  of  inter¬ 
faces  in  motion,  as  compared  to  the  influence  of  the 
hydrodynamic  effects. 

3.  Ihe  developnsnt  of  a  general  theory  of  inmiscible 
fluid  displacement  in  porous  media  idticb  accounts  for 
the  observation  of  viscous  fingering,  (cf •  Rose,  VZZZ) . 

ACaOTOWLBDQMBHT 

The  work  undertaken  and  achieved  in  Contract  CICL-517 
%ras  made  possible  largely  by  the  technical  assistance  given 
by  Messrs.  Buckles,  WUlkow  and  Zeffert  of  the  Chemical  warfare 
Laboratories,  and  by  Professors  A.  E.  Scheidegger,  John  R. 
Philip  and  Zrving  Fatt  during  their  respective  periods  as 
visiting  professors  at  Zllinois. 


SYXnRESZS  OF  THE  DYNMUCS  OF  CAPZLLMY  ZNBZIZTZOV 

By 

Whiter  Rose 

1)  .  Given*  a  single  (peripher^ly  closed)  pore  channel  of 

specified  volume*  V*  defined  by  the  equation 
S  »  S  (x*y*s) 

Thus*  the  space  cooanllnates  of  all  points  lying  on  the 
surface  of  the  pore  channel  are  Xnown;  moreover*  the 
pore  channel  will  have  entry  and  exit  "ends"  as  defined 
by  the  closed  space  curves  tdiich  superimpose  as  terminal 
positions  on  the  surface  8.  (See  figure  1  ); 

2)  .  As  an  initial  condition*  let  us  assume  that  the  pore 

channel  contains  some  volume  of  wetting  fluid*  say 
which  is  smaller  than  V*  and  is  elsewhere  filled  with 
a  non%ietting  fluid  (e.g.  the  saturated  vapor  of  the 
wetting  liquid)  *  so  that  it  can  be  said  (at  all  sub* 
sequent  times  as  well  as  initially) : 

V  •  V  V 

w  n 

To  further  speci:^  the  initial  condition*  it  is  assuawd 
that  both  the  entry  and  exit  ends  of  the  channel  are 
temporarily  closed  by  iasMiginary  barriers*  and  that  a 
sufficiently  long  period  of  time  has  elasped  so  that 
the  fluids  have  asstmed  a  spatial  configuration  which 
makes  the  free  surface  energy  of  the  system  a  minisnim 
value. 

The  first  problem*  therefore*  is  to  specify  this  minimum 
energy  configuration  by  finding  the  equation  of  the 
fluid-fluid  interfacial  surface.  This  solution  in  the 
limit  will  also  give  the  closed  space  curve  which  defines 


or. 
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th«  contact  of  the  fluid-fluid  Interface  with  the 
surface,  8. 

The  free  surfeice  energy  of  the  pore  channel  system 
with  Its  contained  fluids  (for  all  configurations)  1st 
B-X  "•sw,  sn,  wn 

idtere  the  Y '  s  are  the  Interfacial  tensions  (ergs  per 
unit  area) ,  and  the  ^ '  s  are  the  surface  areas  of  the 
two  solid-fluid  and  the  one  fluid-fluid  Interfaces. 

A  second  useful  relationship  Is  Yotmg' s  equation, 
namely: 

%dilch  says  that  at  equilibrium,  the  angle  of  contact, 
made  by  the  fluid -fluid  Interface  with  the  surrounding 
pore  channel  wall  will  have  a  fixed  and  Invariant  value. 

Now,  since  the  ftinctlon,  S,  la  known,  the  total  surface 
area  of  the  pore  cheumel,  0^,  can  be  calculated.  Also 
to  be  noted  Is  the  fact  that 

(fr  "  ^w  ^  ^^n 

Thus,  the  configuration  of  minimum  surface  energy  Is 
fotind  by  applying  the  methods  of  the  calculus  of  variations. 


setting  " 

0,  which  reduces  to: 

"  (r«(r~ 

*  C.  'T* 

-  <ir«. ' 

4r». 

In  the  calculation  being  considered,  gravity  (and  other 
body  fotces)  are  assumed  to  be  absent;  therefore,  a  further 
simplifying  condition  can  be  Imposed,  namely,  that  the 


*sw  designates  the  solld-wettlng  fluid  Interface,  sn  designates 
the  solld-nonwettlng  fluid  interface,  wn  designates  the 
%#ettlng-nonwettlng  fluid  Interface. 
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•qulIibriuB  configuration  of  tha  fluid*>fluid  Interface 
will  be  a  conetant-curvature  evirface,  or: 

{1/r^)  +  (l/r^)  •  K 

where  K  is  the  constant  curvature  of  the  interface  at 
every  point  aiid  the  r’  a  are  the  principal  radii  of 
curvature  at  every  point. 

In  general,  it  nay  be  expected  that  if  no  path  of 
approach  is  postulated,  there  will  be  one  (and  only  one) 
configuration  of  ninimnn  surface  energy  which  exists 
conslstant  with  the  geonetry  of  the  pore  channel  (as 
defined  by  the  function,  S,  and  the  parameter  V) ,  with 
the  initial  value  of  V  ,  and  with  the  constants  of  the 
system  (i.e.  the  ^*s  and  the  value  of  0^)  •  That  is, 
no  matter  how  the  wetting  liquid  (as  measured  by  the 
volume,  V^)  is  introduced,  paths  via  diffusion*  exist 
when  the  time-period  Is  long  \diich  malce  it  possible  for 
a  truly  minimal  energy  condition  to  be  achieved.  (NCte, 
if  the  approach  to  equilibrium  were  by  a  flow  path  alone, 

^  one  could  expect  in  the  general  case  that  energy  "barriers" 


*  and  capllleury  condensation  (eulsorptlon) 


would  b«  oncountorodf  which  is  to  way  that  a  aataatabla 

■inlanin  anargy  configuration  %«ould  ba  attainad  ralatad 

,  * 

to  tha  approach  path)  . 

Iha  final  atap  at  'zatm  ataga  of  tha  calculation  ia  to 
derive  the  curvature  of  the  fluid-fluid  interface  fron 
a  detarsiination  of  tha  aquation  of  tha  wn-intarfaca« 
frcMB  which  ia  then  deduced  tha  aiagnituda  of  tha  praaaura 
diacontinuity  aaaociatad  with  the  interfacial  curvature 
by  the  Laplacian  relationahip: 

A* 

3) .  We  next  aaaune  that  an  infinitaaimal  additional 
anount  of  itfetting  fluid  is  introduced^  aa  aeaaurad  by 
A  V  ,  ao  that  the  new  volune  of  wetting  fluid  ia  (V  )* 

~  W  W 

For  the  purpoaea  of  conceptual  viaualiaation*  it  will  ba 
deairable  to  have  United  in  the  firat  place  conaidaration 
to  that  United  claaa  of  pore  channel  geonetriea  characteriaed 


«  The  analogy  of  a  ball  "rolling  down  a  pocketed  hill"  will 
clarify  thia  point.  The  kinetic  anargy  of  tha  ball  ia  dia- 
regarded,  and  the  potnoitial  energy  of  tha  ball  ia  taken  aa  e 
correaponding  to  the  aurface  energy  of  the  capillary  channel 
ayatan.  Zn  thia  caae,  tha  ball  rolla  to  tha  firat  pocket 
(energy  barrier)  and  atopn>  long  aa  ita  approach  to  "aquilibriun" 
ia  linited  by  aaying  that  tha  approach  path  of  tha  ball  nuat 
be  naintained  on  the  aurface  of  the  hill.  To  bypaaa  the  trap, 
external  work  nuat  be  done;  otherwiaa,  an  independent,  aacondary 
path  of  approach  to  a  lower  equilibriun  nuat  ba  poatulatad. 

Zn  the  capillary  tuba  ayaten,  vapor  diffuaion  and  capillary 
condenaation  afford  thia  indapandant  path.  That  ia,  tha  inter¬ 
face  can  apontaneoualy  nova  forward  by  flow  aa  long  aa  tha  inter¬ 
face  energy  of  the  ayaten  ia  alwaya  dacraaaing.  Tha  flow  than 
stopa  (aay,  tanqporarily)  whanavar  a  further  forward  novenant 
meana  that  tha  anargy  of  tha  ayaten  ia  incraaaad  (i.a.  tha 
ball  ia  in  the  botton  of  a  pocket) ;  however,  if  a\if ficiant  tine 
ia  given,  tha  watting  fluid  can  ba  tranaportatad  ahead  via 
diffuaion  capillary  condenaation,  to  provide  for  a  conplata 
approach  to  a  atata  of  nininun  aurface  anargy.  Tha  forward 
wetting  fluid  nay  or  nay  not  later  coalaace  with  tha  wain  body 
of  wetting  fluid. 


•o  that  the  equillbrlua  configuration  of  tha  initial 
voluna*  of  the  wetting  fluid  ie  located  (at  least  in 
part)  continuous  with  the  so<«alled  entry  end  of  the 
channel. 

In  any  case>  we  now  seek  to  find  the  new  (BiniaROD  sur¬ 
face  energy)  configuration  of  the  new  volvns  of  wetting 
fluid  (i.e.  ^  nethods  outlined  in 

paragraph  two  above.  In  this  case,  however,  it  will  be 
desirable  to  limit  first  attention  to  the  configuration 
for  which  an  approach  path  (in  the  flow  sense)  is  obviously 
extant;  however,  ianediately  a  great  difficulty  is  encountered. 
For  example,  even  though  ^  is  made  very  small 
(approaching  zero),  the  very  act  of  this  addition 
necessarily  pushes  the  fluid -fluid  interface  forward, 
after  *Aiich  a  minimum  surface  energy  configuration  is  to 
be  sought.  However,  in  the  general  case,  the  geometry 
of  the  pore  space  (as  identified  by  the  function,  S) 

Just  ahead  of  the  initial  position  of  the  fluid-fluid 
interface  may  be  such  that  external  work  will  have  to  be 
done  (to  an  underminable  extent)  to  exceed  the  threshhold 
value  of  scxne  incipient  energy  barrier. 

The  foregoing  is  a  problem  since  it  is  the  ultimate 
objective  of  this  work  to  describe  the  spontaneous 
capillary  iadKlbition  of  a  wetting  fluid  into  a  pore 
channel  of  given  geometry.  This  is  to  be  done  by 
noting  the  curvature  of  the  fluid -fluid  interface 
after  each  incremental  addition  of  a  volume  of  wetting 
fluid,  V^,  from  tdiich  can  be  derived  the  time-changing 
value  of  grad  P.  The  latter  is  then  utilised  as  the 
driving  force  term  appearing  in  the  Navier-Stokes  equation, 
which  upon  integration  gives  the  rate-of-interface  movement. 


In  capillary  flows  of  the  type  helng  considered^  gravity 
(and  other  body  forces)*  and  press\ire  forces  sMty  be 
acting  together  with  the  surface  forces  to  produce  notions 
of  interest.  It  seens  necessary  however*  to  first  deduce 
the  aagnitude  and  influence  of  the  capillary  (surface) 
driving  force  as  acting  alone  as  these  relate  In  a  conplex 
way  to  the  function*  8.  Once  they  are  deteznined*  then 
the  added  effect  of  supplying  pressure  energy  and  of 
considering  the  action  of  external  body  forces*  san  be 
accounted  for  by  simple  superposition. 

It  seems  desirable*  therefore*  to  limit  attention  first 
to  a  pore  energy  geometary  such  that  the  function*  8* 
will  define  a  surface  idK>se  cross-sectional  perimeter 
will  Increase  continuously  when  proceeding  from  the  entry 
end  (where  Initially  the  wetting  fluid  wiil  be  located) 
to  the  exit  end*  Indeed*  because  of  the  Above  considerations* 
and  to  achieve  additional  simplicity*  we  may  say  that  it 
will  be  dcwirable  to  first  limit  attention  to  a  capillary 
channel  tdiich  is  symmetrical  around*  say*  the  s-axis* 
and  whose  radial  coordinate*  say  R*  increases  continuously 
from  z"0  to  (where  L  is  the  length  of  the  channel) . 

With  the  above  sis^lification*  the  surface  8  will  be 
characterised  by  the  function*  R  ■■  R(s)*  so  that  when 
dR/ds  is  small*  large  values  of  ^  considered* 

and  when  dR/ds  is  large*  the  values  of  to  be  con¬ 

sidered  will  be  correspondingly  reduced  to  gain  in 
resolution  of  the  tlme*changing  magnitude  of  the  capillary 
driving  force. 
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FLUID-FLUID  I2ITERFACES  IN  UNSTEADY  MOTION 

By 

Walter  Ro0o,  Harry  Fara  and  N.  Chaudharl 

(1; 

In  a  recent  paper  reference  was  made  to  bow  a  con¬ 
stant  curvature  Interface  could  be  expected  as  defining  tbe 
Junction  of  two  limlsclble  fluids  in  steady  motion  through 
a  capillary  tube  of  constant  cross-section.  With  this 

O'!  A 

assumption^  Rose  and  Heins  toolc  the  contact  angle, 
as  a  measure  of  the  curvature  of  advancing  interfaces,  and 
they  reported  data  for  the  steady  flow  of  certain  liquids  in 
circular  capillaries  Which  seemed  to  fit  the  linear  relationship: 

COB  0^-  coa  a; 

where  the  i^idependent  variable,  Z  ,  is  the  velocity  of  inter¬ 
face  movement,  the  slope,  characteristically  is  negative, 
and  the  intercept,  cos  /y  ^  is  the  cosine  of  the  equilibrium 

contact  angle  as  defined  by  Young’ s  equation.  We  now  prepare 

Oi 

to  test  Equation  (1}  against  the  data  of  LeQrand  and  Rense 

Ibe  contact  angle  rate-dependency  given  by  Equation  (1) 

(41 

has  a  font)  conslstant  witli  the  early  conclusions  of  Ablett 
and  the  later  speculations  of  Barter  j  indeed,  Siegel 
indepeculently  has  made  use  of  the  equivalent  of  Equation  (1) 
to  explain  some  of  his  findings.  Since  values  of  cosine  ^ ^ 
less  than  -1  are  not  allowed,  ho%i)ever  Equation  (11  in  the 
niost  favorable  cases  will  only  be  descriptive  of  observed 
data  at  low  values  of  Z" 

It  is  our  purpose  now  to  show  that,  as  a  useful 
approximation.  Equation  (1<  c2U)i  be  expanded  in  tbe  form  of  a 
Maclaurin  Series  when  acceleration  as  well  as  velocity 
deftendency  is  to  be  taken  into  account.  Thus: 


^2} 


cos  =»  cos  +  HZ’  *  HZ-  + 


A 

(h) 


^  J+  t> 


u 


(2) 


vdtore  the  Independent  verlable,  Z*,  la  the  acceleration  of  the 

Interface  noveaient,  and  the  M-co~e£flclenta  are  the  tlae-sero 

Halts  of  the  partial  derivatives  of  cos  0  with  respect 

to  the  associated  Independent  variables.  S'  and  Z*. 

Let  us  consider  vertical  rise  an  Incoapresslble  of 

wetting  liquid  in  circular  capillaries  In  the  aanner  of 
(7^ 

Brittin  (that  is,  neglecting  the  viscosity  and  Inertia 
of  the  vapor  ahead  of  the  iiobiblng  liquid,  neglecting  the 
entry  turbulence,  employing  the  Bolseullle  steady^state  viscous 
resistance  term  in  the  unsteady ‘•^state  systea,  etc.K  We  aay 
then  write: 


(Z  2"  )  ’  +  A  (Z  Z'  ^  +  B  Z  +  C  cos  ^  "  0  (3) 

vribere  A  “ 

B  »  g  (i.e.  the  acceleration  due  to  gravity) 

C  -  -*(2  T  R) 


and  %d)ere  Is  the  viscosity,  ^  is  the  density,  and  T  is 
the  surface  tension. 

Talcing  only  the  first  three  terms  of  Equation  (2)  as  an 
approximation  and  conblnlng  with  Equation  (3) ,  gives  an 
integration: 


(Z  Z’  )  +  (A/2)  (z  )  4-  B  /z  dt  « 


Jz  dt  -  c  ^ 


cos  ^  ,t  +  M  Z  +  M 

B  V 


With  the  aid  of  Equation  4,  then,  the  constants,  and  M^, 
can  be  evaluated  from  a  set  of  experiaental  data  of  Z  as  a 
function  of  time,  t,  by  a  leaat>squares  analysis  Involving 
the  simultaneous  solution  of  the  following  set  of  cKpiations: 


(4) 


(5) 


'3} 

«  •  ^  +  "v  +  ®iV 

®  *  -V  +  »iV  <»i> 

whttre  Is  the  value  of  the  left-hand  aeinber  of  Equation  (4) 

ait  tine  t^,  and  la  the  ratio  of  Z'  to  Z  at  tine  t^. 

In  Table  I  we  give  the  values  of  N  and  M  «diich  have 

V  a 

been  calculated  by  Equations  (5)  fron  the  data  obtained  la 

the  nine  capillary  rise  experinents  ireported  by  LeOraod  and 
(3)1 

Rense  .  Ihe  following  points  are  of  interest: 

1)  .  Values  of  characteristically  are  one-tenth  the 

corresponding  values  of  therefore*  one  apparently 
is  Justified  in  dropping  the  fourth  and  higher 
terms  of  Equation  (2) . 

2)  .  In  a  decelerating  system*  the  contact  angle  is 

greater  at  a  particular  velocity  than  if  the  flow 
were  a  steady  flow  at  the  same  velocity,  therefore* 

M  is  characteristically  negative  auod  character¬ 
istically  is  positive*. 

3) .  During  the  initial  ia^lslvc  motion  as  Z*  goes  from 

zero  to  some  maximum  value*  Z”  is  positive*  so  that 
the  term*  M^Z*  serves  as  a  check  on  the  term*  M^Z*  * 
and  unreasonable  valties  of  cosine  are  not  predicted 
by  Equation  (2).  N.B.  Compensation  also  afforded  by 
other  neglected  terms. 

4)  .  seems  to  correlate  with  the  dimensionless  group 

formed  as  the  ratio  of  the  viscous  to  capillary 
forces.  That  is*  the  greater  the  viscosity  (other 
things  being  equal) *  the  greater  will  be  and* 
conversely*  the  greater  the  surface  tension  (other 
things  being  equal)  the  smaller  will  be  N^. 

*  Indeed*  it  is  as  if  the  decelerating  interface  "remembers” 
the  higher  previous  value  of  contact  angle. 


f4J 

This  note  on  uneteedy  aotlone,  like  the  previous  one 
(IS 

on  steady  atotions  ,  in  a  heuristic  way,  gives  infomation 
about  the  curvature  of  fluid-fluid  interfaces  amving  in 
capillary  channels. 

TABLE  I 

LeGrand  and  Rense 
Experiment  Number 


3  (water  in  cm, 

tube)  -  O0O37  Mg  >•  0.0020 


4  (EtOH  in  0-0242  «m 

tube)  *  0.028  M,  -  0,0018 


The  complete  table  is  not  presented  at  this  time  inasmuch 
as  the  values  obtained  by  desk  calculator  operating  on 
graphically  determined  values  of  Z”  and  on  plamimeter-values 
of  the  integral  of  Zdtp  were  not  showing  the  desired  resolution, 
Inddue  time;  a  machine  programmed  solution  of  Equation  5  will 
be  available.  To  be  noted;  none  the  lessp  is  the  fact  that 
the  values  given  above  correspond  closely  to  thos»  given  in 
Tables  I  and  II  of  the  13th  Supplementary  Statement  of  this 
Progress  Reports 
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VISCOUS  FINOERmS  AS  A  FEATURE 
OF  STABLE  DISPLACENEMT 

By 

Walter  Rooe 

1  Historical  Conaiderations 

2  Microscopic  Dlaplacement  Theory 

3  Problems  in  Macroscopic  Representation 

4  Oausl-lnstablllty  analogies 


Abstract 


In  recent  years  it  has  become  popular  to  suppose  that 
the  stability  of  a  displacement  fro^\t  is  determined  by  the 
magnitude  of  the  so-^called  mobility  ratio  and  other  associated 
factors.  While  agreement  has  not  been  obtained  on  all  points 
the  summaries  presented  by  Scheldegger  (19601  and  Collins (1961) 
may  be  taken  as  the  representative  views  which  current  book 
authors  choose  to  advance. 

In  this  Note,  we  take  exception  to  these  popular 
theories  «Ailch  account  for  the  observation  of  viscous  fingering 
as  a  manifestation  of  Instability  phenomena.  First,  we 
emphasize  that  it  is  a  viscosity  ratio  condition,  and  not  a 
mobility  ratio  condition,  which  must  be  met  before  viscous 
fingering  will  be  observed.  Second,  we  emphasize  that  the 
flood  front  Is  not  a  well-defined  surface*  to  vdilch  Instability 
considerations  are  to  be  applied,  but  rather  It  Is  simply  a 
moving-zone  region  having  a  high  density  of  microscopically 
distributed  fluid-fluid  Interfaces,  which  region  is  marked  by 
a  high  normal  saturation  gradient.^  Third,  we  emphasize  that 
the  physico-chemical  Interactions  and  the  laws  of  motion  which 
apply  to  the  description  of  displacement  systems^,  at  least  In 
the  microscopic  sense,  are  no  different  for  high  viscosity 
ratio  (and/or  mobility  ratio.:  systems  than  for  low  viscosity 
ratio  (and/or  mobility  ratio)  systems,  therefore,  the  analysis 

^"iNor"°«cajnple,  analogous  to  the  Hele-Shaw  cell  interface 
discussed  by  Saffroeut  and  Taylor  (1958) . 


and  dMcriptlon  of  vlacous  flxigerlng  antaila  no  spaclal 
difficulties  beyond  those  net  in  ordinary  (so-called  stable) 
displacenmt  problems. 

As  for  tile  first  point,  eoqperinent  shows  that  viscous 
fingering  is  never  Observed  unless  the  viscosity  of  the  entering 
fluid  is  leas  than  that  of  the  replaced  fluid,  no  natter  idiat 
value  the  nobility  ratio  nay  have  had  during  the  experinent. 

As  for  the  second  point,  it  can  be  shown  that  the  wave  length 
of  incipient  fingers  corresponds  to  the  wave  length  of  nicro- 
heterogeniety,  and  both  involve  linear  dinensions  coinparable 
to  those  idiich  neaeure  the  areial  extent  of  the  various 
fluid-fluid  interfaces  in  notion  (i.e.  a  few  pore  dianeters 
at  the  nost  between  contiguous  pore  wall  contacts) ;  therefore, 
once  can  talk  only  about  a  condition  of  quasi-instability  in 
analogy,  for  exanqple,  to  the  unstable  configuration  of  the 
Interface  seen  in  the  Hele-Shaw  cell,  and  in  niscible  fluid 
porous  media  displacement  systens  %dien  the  entering  fluid  is 
less  viscous  than  the  replaced  fluid.  Finally,  as  for  the 
third  point  above,  we  show  herein  that  the  Buckley-Leverett 
procedure  for  describing  the  unsteady  states  of  lamiscible 
fluid  nixture  flow,  can  be  applied  to  represent  viscous  fingering 
with  the  same  degree  of  success  as  experience  shows  it  applies 
to  the  so-called  stable  displacement  processes. 
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Physics  of  Inolsclble  Fluld*Fluid  nisplacsnent 

By 

Walter  Rose 


1.  Introduction 

Hie  physics  of  flow  through  porous  itedla  is  a  subject 
of  tecognlsed  Inportemce  and  wide  popularity.  It  Is  also  a 
subject  of  great  complexity,  and  even  the  recent  book 
(Polubarlnova~Kochlna,  1952;  Carman,  1956;  Boupeurt,  1957; 

Lykow,  1958:  Scheidegger,  1960;  Collins,  1961J  acknowledge 
that  there  Is  still  much  to  be  learned.  In  any  case, 
view  is  adopted  here. 

To  visualize  the  problem  fully,  one  must  become  acquainted 
with  the  general  features  of  transport  phenomena  aa  discussed 
for  example  by  Bird  et  al  (.I960) ,  and  with  the  particulars  of 
tlieoretlcal  fluid  mechanics  as  treated  for  example  by  Landau 
(1959> .  With  tills  foundation,  one  can  begin  to  construct  the 
equations  of  motion  whose  integrals  are  wanted  to  predict  the 
details  of  porous  media  flows. 

It  is  the  purpose  of  this  paper  to  show  that  even  idten  we 
limit  attention  to  rather  simple  and  elementary  cases,  great 
difficulties  in  analysis  often  are  encountered.  These  may 
occur  because  the  porous  media  geometries,  in  general,  are  not 
specifiable;  moreover,  the  large  characteristic  surface  area 
of  the  pore  structure  gives  rise  to  unusual  viscous  drag  effects 
and  to  unique  capillary  (fluid •fluid-solid)  interactions.  In 
consequence,  there  is  a  great  temptation  and  tendency  to  over-" 
simplify  the  mathematical  model  chosen  to  represent  the  porous 
media  flow  system,  leading  to  the  categories  of  uncertainty 
already  dlscussod  with  great  insight  by  Masse  (1958). 

Let  us  limit  attention  arbitrarily  to  some  sioqple  diS'< 
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placanent  problem  t«here  one  fluid  is  seen  to  replsce  another 
immiscible  fluid  from  the  pore  space  of  an  interconnected  net- 
mork  of  capillary  channels  (i.e.  a  porous  medium).  The  widely 
accepted  solution  to  this  problem  is  that  given  in  ^  2  of  this 
paper«  which  solution,  however,  is  believed  to  be  suspect  on 
three  counts,  namely: 

1.  No  explicit  account  is  taken  in  the  general  theory 
as  referring  to  the  onset  of  instabilities  believed  by 
some  to  develop  at  the  later  stages  of  displacement  (i.e. 
when  the  mobility  of  the  replacing  fluid  finally  becomes 
higher  than  that  of  the  replaced  fluid) .  Ihls  matter  is 
an.ly*«d  in  ^3  of  this  paper. 

2.  n:e  general  theory  is  derived  from  the  untested 
assumption  that  the  unsteady-states  of  multiphase  flow 
can  be  described  by  reference  to  a  Darcy  law  analog.  In 
^4  this  assumption  is  analysed  and  then  abandoned. 

3.  The  general  theory  treats  the  action  of  capillary 
forces  in  a  way  «ihlch  is  unnecessarily  inexact,  as  dis¬ 
cussed  in  §S  of  this  paper. 

The  three  topics,  enumerated  above,  axm  introduced  here 
heurlstically  inasmuch  as  they  have  been  largely  Ignored  by 
previous  authors.  For  example,  there  is  no  mention  of  them  in 
the  books  cited  above.  While  it  must  be  admitted  that  the 
general  theory  which  ignores  these  points  gives  good  results 
in  some  cases,  it  will  be  shown  that  modifications  must  be  in¬ 
troduced  to  predict  correctly  the  course  of  other  practical 
cases  of  porous  media  flow. 
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1.  lk>  •xpliclt  account  is  takan  in  tha  ganaral  thaory  which 
will  rapraaant  tha  onaat  and  growth  of  viscoua  flngaring  which 
is  saan  to  occur  %dien  certain  conditions  are  aet  (including 
the  condition  that  the  replacing  fluid  is  nore  mobile  than  tha 
displaced  fluid.  It  is  true  that  guasi-insthbility  considerations 
have  been  introduced  by  soma«  but  in  ^  3  of  this  paper  it  is 
shown  how  the  latter  approach  begs  rather  than  resolves  the 
question. 

2.  Statement  of  the  Problem 

Suppose  we  have  a  porous  medium  with  given  boundary 
conditions  imposed,  and  with  given  initial  conditions  of 
saturation  and  saturation  distribution  of  two  (or  more) 
imniscible  fluid  saturants.  The  conditions  are  chosen  so 
that  the  energy  and  mass  transport  affects  are  unisqportant, 
therefore,  we  can  limit  attention  to  momentum  transport  phenomena. 
In  particular,  we  seek  a  way  for  describing  (predicting)  the 
course  of  events  during  the  ensuring  unsteady-states  marked  by 
saturation  changes,  as  a  function  of  time  and  position, 
throughout  the  porous  medium  system. 

Without  serious  loss  of  generality,  tha  precise  statement 
of  the  problem  can  be  greatly  simplified  at  the  onset,  as 
follows : 

1.  By  choosing  an  Isotropic  porous  maditun  in  tha  sense 
that  tha  porosity  and  permeability  (l.a.  the  interstitial 
geometry)  are  homogeneously  uniform  throughout  the  space 
of  interest  (at  least  in  all  flow  directions),  and  by 
saying  that  the  internal  pore  surface  is  continuous  and 
chemically  the  same  at  every  point. 

2.  By  stipulating  that  the  solid  matrix  itself  is  non- 
porous,  rigid.  Incompressible,  insoluble  and  inert. 
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3.  By  having  only  tvo  fluids  as  the  pore  saturants« 
one  of  which  is  taken  as  the  >wtting  liquid  (that  Lm, 

the  equilibrium  contact  angle  measured  through  this  liquid 
is  less  than  90" ) ,  and  by  saying  these  imsiisclble  fluids 
are  Newtonian,  each  being  a  single  component  phase 
mutually  satxurated  with  the  other,  and  each  having 
negligible  compressibility  but  having  comparable  (say, 
equal)  viscosity  and  density. 

4.  By  limiting  attention  in  the  macroscopic  sense  to 
one>-dimensional  flow  (e.g.  the  source  and  sink  are  the 
opposite  faces  of  an  otherwise  encased  cube) ,  udiere  the 
Reynolds  Nuadder  is  low  enough  (say,  leas  than  unity) 

so  that  turbulence  is  avoided  and  so  that  the  nonlinearities 
in  streamline  flow  will  be  minimised. 

Obviously,  the  system  under  consideration  could  be  con¬ 
structed  by  taking  a  macroscopic  cubic  element  of  a  random 
packing  of  (say,  well-sorted  glass)  spheres  as  the  porous 
medium,  emd  by  taking  (say)  pure  water  and  some  selected  hydro¬ 
carbon  fraction  as  the  fluid  pair.  As  a  simple  condition,  we 
could  let  the  medium  initially  be  saturated  completely  with 
the  wetting  liquid,  so  that  our  porblem  is  to  trace  the  history 
of  %iihat  happens  when  the  nonwetting  fluid  is  forced  to  displece 
and  replace  the  wetting  liquid  (drainage) ,  and  subsequently  to 
trace  the  history  of  udiat  happens  when  the  wetting  liquid  is 
brought  aack  into  the  porous  medium  (imbibition) .  n»e  latter, 
of  course,  can  occur  spontaneously  (via  the  action  of  capillary 
forces,  with  or  without  additional  press  ^re  forces  externally 
iaqposed) .  In  any  case,  experience  showid  ^  that  repeated  cycles 
of  draine^r*  end  imbibition  displacement  give  rise  to  hysteresis 
effects  of  interest,  and  we  shall  want  our  fheoary  to  account 


*cf.  scheidegger,  loc.  cit. 
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for  this  irrsvsrslbillty  feature. 

Much  of  the  pertinent  work  on  the  problen  at  hand  is 
cited  by  Scheidegger  (1960) .  The  first  attesipted  quantitative 
statenent  'was  given  by  Richards  (1931) ,  but  siany  other  in¬ 
vestigations  have  been  isiportant  in  the  evolution  of  the  con¬ 
cepts  which  now  are  to  be  exanined  (cf .  Buckley  a.vl  Leverett, 
1942;  Brinkman,  1948;  Klute,  1952;  Philip,  1957;  Cardwell,  1959; 
Blair  ^  al.  1958;  Bovanessian  and  Payers,  196  ) . 

The  analyses  necessarily  starts  with  a  mass  continuity 
statement  of  the  type: 

(2.1)  div  (  vj  -  -  P^  (S,  pj 


'i  f  V 


from  which  can  be  derived  (Blair  et  al.  1958) : 


(2.2) 


/•  L  3/^ 


ICfy  ^ 


''7 


on  assuming  that  Darcy*  s  law  applies,  and  that  the  effective 
permeability  and  the  capillary  pressure  are  known  functions 
(e.g.  empirically)  of  saturation  according  to: 

(2.3)  k^  »  k^  (Sj^3  -  -  v^ 

grad 

and 


(2.4) 


P  =•  P  {8  )  P  *  p  p 

Equation  (2.2),  vihen  condoined  with  (2.3)  and  (2.4), 
implicitly  has  the  wanted  solution: 

(2.5)  (x,t) 

and  thus.  In  the  common  view  (cf.  Hovanesslan  and  Payers, 
(19^/  ),  It  represents  "....the  fundamental  differential 

equation. ../«  ^whose^  . solution  gives  a  complete 

description  of  the  physical  problem  being  studied....* 

In  the  statement  of  (2.2)  to  express  linear  displacement 
In  the  absence  of  gravity,  however,  the  following  conditions 
and  conventions  have  been  Imposed,  namely: 

A  V  +  V  »  const. 

/  w  n 


P  *  P  ^  P 

*^n  w 

S  +  S  «  1 

w  n 


(2.6)  /  p^/ y  ■  » 

»  const.  axAJLfJ  »■  0 


Thus,  it  Is  seen  that  the  solution  (2.S)  Is  had  when  the 
saturation  Is  known  as  a  function  of  time  and  position  through-¬ 
out  the  system,  which  means  that,  implicitly,  the  other  depen¬ 
dent  variables  (v.,k  ,p  ,p  )  also  are  determinable  and  known 

X  X  jL  C 

In  terms  of  x  and  t,  while  all  other  parameters  are  given  con¬ 
stants  of  the  system. 

In  the  most  general  case  of  three-phase  saturated  media, 
we  have  In  addition  to  the  four  Independent  variables  (x,y,z  ax^  t) , 
18  dependent  variables  (namely,  three  each  of 
and  three  caplllajry  pressure  parameters  for  the  L2,  13  xmd 
2-3  fluid  parrlngs) ;  therefore,  %#e  need  a  set  of  18  simultaneous 
equations  each  of  the  continuity  statement,  the  equation  of 


•tatei  Darcy's  law.  the  K^<S^)  function;  and  the  capillary 
pressure  definition;  two  equations  cflvlng  p  (s  } ;  and  one 

C  jL 

equation  of  the  form  ^  *  1) . 

Indeed^  claims  have  been  made  (  Uou^las  £t  ai,  1959  3 

that  more  general  solutions  can  be  obtained  of  the  form: 

(2.7)  <x,y,z,t) 

and,  to  the  extent  that  machine  programming  and  the  supple¬ 
mentary  exp^erlmental  work  (to  give  the  functions  2.3  and  2.4) 
can  be  undertaken  with  reasonable  labor,  ultimately  It  may  be 
possible  to  treat  multl-dlmenslonal  flows  with  three  or  more 
(miscible  and/or  Immiscible)  fluids  present*,  taking  gravity 
as  well  as  capillarity  Into  account,  for  transport  when  con¬ 
ditions  opposite  to  (2.6)  apply.  In  porous  media  which  are 
neither  Isotropic  nor  homogeneous.  Consideration  of  these 
multiple  possibilities,  however,  adds  nothing  to  our  present 
analysis,  except  thereby  we  shall  draw  the  Inference  that  It 
is  self --defeating  to  increase  the  order  of  complexity  of  a 
problem  until  the  elementary  solutions  have  been  fully  evaluated. 

In  point  of  fact.  It  is  our  Immediate  purpose  to  see  In 
(2.2)  is  a  proper  representation  of  the  dynamics  of  displacement, 
that  us,  tdiethar  or  not  one  can  avoid  the  Indetermlnancy  of 
a  partial  differential  equation  with  three  dependent  variables 
(l.e.  2.1)  by  reducing  the  number  to  one  (l.e.  by  operating  on 
2.2  with  2.3  and  2.4).  It  Is  only  In  this  way  that  the  solution 
(2.5)  can  be  obtained,  therefore,  the  availability  and  accessablllty 
of  solutions  of  the  type  (2.7)  Is  a  contingent  matter. 

It  Is  necessary.  In  consequence,  to  examine  In  detail  the 
meaning  and  applicability  of  Darcy's  law  (2.3)  and  the  capillary 
pressure  function  (2.4);  for.  If  these  Intermediate  steps  are 

*In  miscible  fluid  displacement  systems,  a  concentration  term 
Is  introduced  for  the  saturation  term  of  (2.1),  and  a  dlffusivlty 
pareuneter  replaces  the  poroslty-'permeablllty  function  applicable 

to  Immiscible  fluid  displacement. 


rejected  In  part  or  in  «ihole,  the  hope  of  obtaining  even 
approximate  solutions  of  the  type  (2.5)  and  (2.7)  is  corres* 
pondingly  diminished.  Firsts  however,  it  is  important  to  con¬ 
sider  the  quasi-instability  matter  udiich  gives  rise  to  the 
development  and  growth  of  so-called  "viscous  fingers"  ^en  the 
displacing  fluid  has  a  higher  mobili^  than  the  replaced  fluid. 
The  latter  feature  often  is  treated  as  a  separate  question 
(Scheidegger  and  Johnson,  1960) ,  but  it  is  evident  that  viscous 
fingering  is  an  important  aspect  of  the  unsteady-states  of 
miscible  as  well  as  immiscible  fluid  displacement  processes. 

Our  phenomenological  theory  immediately  follo%#8. 

3.  Insthbilitv  Considerations 

If  instabilities  occur,  the  continuity  implied  by  (2.1) 
still  applies;  therefore,  we  must  examine  how  the  instabilities 
will  affect  the  form  of  the  solution  (2.5). 

Let  us  suppose  that  Darcy's  law  (2.3)  describes  the  flow 
transport  of  each  fluid  pahse,  even  though  there  may  be  a  high 
viscosity  contrast  between  leases,  and  even  when  saturations 
are  changing  rapidly  with  time.  We  take  this  latter  as  a 
siii^lifylng  assumption,  tentatively,  until  the  matter  is  dis¬ 
cussed  in  greater  detail  in  ^  4  to  follow.  In  this  %«ay,  we 
are  free  to  discuss  the  criteria  which  give  rise  to  the  onset 
and  growth  of  the  viscous  fingers  (idilch  are  the  manifestation 
of  the  instability) ,  and  to  employ  the  convenience  of  the  02urcy 
law  concept  in  formulating  the  problem  heuristically. 

We  introduce  first  the  mobility  proportionality  factor 
as  a  measure  of  the  ease  with  %diich  a  fluid  phase  can  be  forced 
to  move  through  a  porous  medium.  Obviously,  fluid  mobility 
is  a  function  of  the  saturation  and  the  saturation  distribution 
configuration.  In  the  Darcy  law  frame  of  reference,  the  nobility 


(9) 


can  be  expressed  as  an  effective  permeability  to  viscosity 
quotient,  dependent  on  tine  and  on  the  space  co¬ 

ordinates,  as  well  as  on  the  past-history  (hysteresis  of  satu¬ 
ration  change)  of  the  systen  in  a  way  idilch  corresponds  to  the 
meaning  of  A  mobility  ratio  therefore  can  be  defined 

arbitrarily  as: 

(3.1)  M  -  ^2^ 

%«here  the  convention  Is  adopted  that  fluid  1  Is  the  displacing 
(entering)  fluid,  and  fluid  2  Is  the  displaced  (replaced) 
fluid. 

We  now  Introduce  the  hypothesis  that  for  horlsontal  flow 
(l.e.  no  gravity  forces  considered  to  be  acting),  the  moving 
flood-front  will  be  unstable  for  all  sensibly  finite  forward 
velocities  when  M  Is  sensibly  less  than  \mxty.  This  postulate 
can  be  justified  qualitatively  by  considering  the  instability 
criteria  of  Chuoke  et  aX  (1958) ,  and  by  noting  that  tie  critical 
wave-lengtii  of  flood-front  pertubatlons  «hlch  must  be  exceeded 
If  fingering  is  to  ensue  numerically  Is  indicated  to  be  sensibly 
smaller  than  the  wave-length  Index  of  heterogenlety  (as  defined 
by  Scheldegger,  1960) .  In  fact  the  latter  apparently  Is  true 
even  for  what  otherwise  would  be  considered  as  homogeneous 
media,  as  will  be  shown. 

Following  Smith  fit  al  (1933),  %fho  introduced  the  useful 
concept  that  a  randma  packing  of  spheres  of  uniform  sise  (i.e. 
a  sensibly  homogeneous  medlism)  Mtually  has  the  form  of  inter¬ 
posed  asserbleges  of  blocks  ahere  a  particular  mode  of  packing 
is  in  predominance  (e.g.  cubic  assenhleges  contiguous  to 
orthorhosblc  and/or  rhombdohedral  assenbleges,  etc.),  we  sus¬ 
pect  that  the  lower  limit  for  the  Scheldegger  wave-length  in¬ 
dex  of  heterogenlety  will  be  measured  by  (say)  ten  to  100  grain 


diaatttcrs.  With  it  thartby  likely  that  incipiant  partuxbations 
ara  a  faatura  of  avan  alanantary  honoganaoiis  aiadia,  tharafora« 
as  wa  conaidar  dlaplacanant  in  leaa  boaioganaoua  (aay«  hataroganaoua) 
inadia«  for  axaaqpla,  vAiara  tha  aorting  of  particla  aixa  and  tha 
uniforaity  of  particla  ahapa  haa  haan  dlaiiniahad«  and/or  a 
Biicro-badding  trand  haa  haan  iaqpoaad,  than  tha  Chouka  gt 
criteria  for  inatability  will  ba  met  for  noat  practical  caaaa 
of  interaat,  aa  long  aa  M  ia  laaa  than  unity. 

wa  adopt  tharafora*  aa  approximately  correct,  tha  propo¬ 
sition  that  M  ^1  iaqplias  that  inatabilitias  will  occur  during 
the  diaplacaoiant  process  under  consideration,  being  prepared 
to  abandon  this  qualitative  stataaiant  later  whanavar  tha  critical 
value  of  N  is  wanted  (somewhere  between  unity  and  aero)  which 
quantitatively  marks  the  transition  batwaan  stable  and  unathhla 
flood-front  character. 

To  facilitate  discussion  and  numerical  analysis,  alb<.it 
without  a  loss  in  generality,  lot  us  introduce  for  (2.3): 

(3.2)  -  (S^  -  /  (i  - 


and: 

(3.3) 


kj^  -  (1  -  -  S^)  /  (1  -  -  8^) 


as  approximations  to  %diat  experiment  gives  as  tha  usual  form 
of  the  relative  parmaab&lity  function  (cf .  hose,  1949) .  Than 
M  may  be  raprasantad  by: 

(3.4)  N  - 

yielding  tha  following  critical  values  of  8^(  where  11*1}  and 
where  tha  onset  of  viscous  fingering  is  assuawd  to  occur,  namely: 

(3.5)  8  - 

C 

whan  ^5^  3)  t  ^  0  (i.e.  a  drainage  process) 


and: 


I 


(11) 


(3.6)  S  « 

c 

«li«n  t  0  (i.«.  a&  iaibD»ltlon  procMs). 

th*  laportaace  of  (3.5)  and  (3.6)  Xlos  in  tha  aasuaqptlon 
that  a  drainage  process  and  an  labibition 

process,  corresponds  to  M  ^  X  (that  is,  instabilities  are  to 
be  expected) . 

Now,  it  has  been  Inferred  already  that  (2.5)  can  be 
obtained  from  (2.2)  through  substitution  of  the  eaqpirical 
functions,  (2.3)  and  (2,4),  only  when  the  flood-front  is  stable 
(e.g.  M  ^  1) .  It  is  laqportant  to  develop  the  rationale  for 
this  restriction.  The  continuity  expression  (2.1)  is  always 
valid  (i.e.  "natter  can  neither  be  created  nor  destroyed*)  quite 
apart  fzon  any  instability  consideration;  IDcewise,  we  are 
prepared  at  the  nonent  not  to  quUd^le  idiether  or  not  use  can 
be  nade  of  the  relative  penaeabllity  function  (2.3)  and  the 
capillary  pressure  function  (2.4).  Ihe  point  is  that  as  long 
as  the  flood-front  interface  is  stable,  the  nicro-faeterogeneities 
of  the  systen  can  be  ignored;  and,  conversely,  when  the  inter¬ 
face  is  unstable  the  nicro-heterogeneitiea  suffice  to  induce 
the  pertubations  «d)ich  give  rise  to  the  ^>servatioa  of  the  viscous 
fingering. 

Figure  1  shows  the  variation  of  M  with  as  obtained 
from  (3.4)  for  different  values  of  the  viscosity  ration,  and 
for  the  Indicated  sets  of  assuBq>tions  about  8^  and  8^  chosen 
to  correspond  to  the  hysteresis  in  the  relative  pemeability 
curves  when  lablbition  and  drainage  displaceaMuit  processes  are 
ccMopared.  Ihe  shift  and  relative  values  of  8^  as  defined  Iqr 
(3.5)  and  (3.6)  are  of  IssBediate  interest. 

Let  us  consider  first  an  Isblbition  process  with  the  initial 
conditions  implied  by  Figure  2.  If  the  wetting  fluid  viscosity 
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!•  considerably  greater  than  that  of  the  replaced  fluid,  8 
will  be  close  in  value  to  (1-S  ) ;  therefore^  taking  x,  as  the 

Jvn  X 

level  of  observation,  S  will  be  seen  to  continuously  increase 

at  as  the  flood-front  passes  under  displacesMnt  conditions 

«A)ere  M  is  favorable  for  nost  of  the  tine.  Indeed,  no  instability 

can  occur  until  8  exceeds  8  ,  but  tine-wise  this  does  not 

w 

happen  until  oost  of  the  nonwetting  fluid  has  been  displaced 

fron  the  level  of  x^,  which  is  to  say  that  the  displacenent 

efficiency  of  the  process  will  be  negligibly  affected  by  any 

viscous  fingering  which  nay  ensue  due  to  the  fact  that  k|^  is 

approaching  sone  vanishing  value. 

On  the  other  hand,  for  inhibition  displacenent  processes 

where  8  is  close  to  8  ,  the  onset  of  Instabilities  at  the 

c 

flood -front  will  be  observed  at  the  level  x^  almost  innediately, 
and  a  corresponding  reduction  in  the  displacenent  efficiency 
of  the  process  can  be  expected. 

nie  final  cases  to  be  considered  have  to  do  with  drainage 
displacement  (i.e.  «diere  the  nonwetting  fluid  replaces  the 
«#etting  fluid) ,  and  it  will  be  noticed  again  fron  Pigxure  1  that 
the  proximity  of  8  to  8  (favorable),  or  to  (1-8  ) 

fVf 

(unfavorable) ,  is  tne  factor  which  determines  at  what  stage 
of  the  displacement  process  the  instabilities  are  manifested. 

Ihus  it  is  seen  that  ^dmn  the  displacing  fluid  is  con¬ 
siderably  more  viscous  than  the  replaced  fluid,  fingering  is 
deferrtMl  until  the  last  stages  of  displacemsnt;  therefore,  the 
instability  phenomena  will  not  be  noticed  beyond  (possibly) 
the  observation  of  a  slightly  higher  "irreducible”  saturation 
of  the  replaced  fluid  than  that  limiting  value  attained  after 
extremely  long  periods  of  extremely  slow  rate-of-flood-front 
advance.  It  is  problems  of  this  sort  for  which  aquation  (2.2) 
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hu  been  developed  by  prevloua  «iorkera  aa  the  appropriate 
deacriptlon^  and  for  %dilch  aolutlona  of  the  type  (2.5)  have 
been  obtained  by  ataking  uae  of  the  eaqpirically  daterained 
functlona  (2.3)  and  (2.4). 

4. 

Zn  lieu  of  the  preaentation  of  the  intended  raeMining 
part  of  Section  3  above«  and  of  Sectiona  4  and  5  referred  to 
in  the  introduction  to  thia  paper*,  an  outline  of  the  arguaMnt 
ia  auanariaed  below. 

The  general  theory  of  iaaiiacible  fluid-fluid  diaplace- 

aMnt  aniat  account  for  the  growth  and  developaMnt  of  viacoua 

fingering,  and  it  ia  likely  that  appeal  can  be  ande  to  in- 

atability  conaiderationa  to  obtain  a  phenoaienolical  deacription. 

In  the  aiore  baaic  aenae,  however,  a  aingle  unifying  repreaentation 

nuat  exiat  which  ia  g«ieral  for  the  caae  of  an  "unfavorable*’ 

aa  i«ell  aa  a  "favorable"  viacoaity  ratio.  Ihia  paper  aeeka 

to  find  that  unifying  repreaentation,  by  conaidering  the  dynaaiica 

of  diaplacenent  aa  viewed  with  reapect  to  the  nicroacopiea  of 

the  ayaten.  The  latter  ia  nade  neceaaary,  in  addition,  becauae 

of  the  conaiderationa  previoualy  brought  forth  by  Roae  and 

Channapragada  (cf.  alao  Ilaaa4,  loc.  cit.)  dealing^  with  the 

uncertaintiea  iinplicit  in  the  uae  of  Darcy'  a  law  to  deacribe 

2 

the  unateady-atatea  of  aixture  flow,  and  ,  with  the  fact  that 
the  capillary  driving  forcea  have  their  aource  in  aiicroacopically 
diatribttted  inter facea. 


*to  appear  later 


(14) 


1.  Bird,  R.  B.  jtt  Phanoiiiana.  Wiley  Book,  I960 

2.  Blair«  P.  M.  Trans.  AZME,  v.  213,  p.  96,  1958 

3.  Buckley,  S.E.  and  Leverett,  M.C. ,  Trans.  AZMB,  v.  146, 
p.  117,  1942 

4.  Brinkean,  H.C.,  Applied  Scientific  Research,  v.  Al,  p.  333, 

1948 

5.  Card%#ell,  W.  T. ,  Trans.  AINE,  v.  216,  p.  271,  1959 

6.  Carman,  P.  C.,  Flow  of  Gases  Through  Porous  Media. 

Academic  Press  Book,  1956 

7.  Collins,  R.  £.,  Plow  of  Fluids  Through  Porous  Materials, 

Reixihold  Book,  1961 

8.  Chouke,  e^  al.AIME  paper  1141,  1958 

9.  Houpeurt,  A. ,  Elements  de  Mechnieue  des  Fluides  dans  lea 
Milieux  Poreux.  Paris,  1957 

10.  Bovanessian  and  Payers,  AZMB* J.Soc .Pet.  Bngr.,  ITaroh  1961,  p.32 

11.  Klute,  A.  Soil  Science,  v.  72,  p.  105,  1952 

12.  Landau,  L.O.,  Fluid  Mechanics.  Addison-wesley  Book,  1959 

13.  Lucien  Masse,  The  Displacement  of  Fluids. 

Oniversity  of  OklalxMoa  Sysfposium  on  Fluid  Flow  Through 
Porous  Media,  Borman,  Oklahoma 

14.  Philip,  J.R.,  Soil  Science,  v.  84,  p.  329,  1957 

15.  Polubarinova,P.ya. -Bochina,  The  Theory  of  ground  Water 
Movement.  Moskow,  1952,  English  Translati<m  in  the  press, 

1961,  hy  R.J.M.  Dewiest,  Princeton 

16.  Richards,  L.A.  Physics,  v.  1,  p.  318,  1931 

17.  Rose,  W.  Trans.  AZMB,  v.  186,  p.  Ill,  1949 

Rose,  W.  and  Channapragada,  R.  ,  AZMB  ,  paper  1549- 
presented  Dnver  meeting,  October,  1960. 


18 


(15) 


Sclio Siw  now  Through  PorouQ  t*ff4Aa> 

MacMillan  Book.  2nd  Edition,  1960 

20.  Scheidegger,  loc.  clt. 

21.  Snlth,  W.O.,  at  al.  Physics,  v.  3,  1933 

22*  J.  J)ou^,la8  at  al,  "A  I,''ethod  for  Calculating  liulti^Dirnenaio.ial 
Iraiaiscibl  "Srs^ac -laent",  Transo  AIIW,  v®  216,  po  297  (1959)© 


£jL|th 


ETETtrORKS  REPRESENTZTO  POROUS  MEDIA  PLOH  SYSTEMS 

By 

Walter  Itose 

Abstract* 

We  limit  attention  to  cases  of  streamline  flow  of 
essentially  incompressible  fluid  (s])  being  transported 
through  porous  media  systems  which  are  characterized  by 
varying  degrees  of  inhomogeniety  and  anisotropy^  however^ 
both  steady  and  unsteady  flow  situations  are  to  be  considered. 

As  a  foundation,  tie  then  state  (in  a  brief  but  rigoroixs 
fashion)  the  consequences  of  applying  Mavier-Btokes  theory 
to  gloe  a  description,  in  microscopic  detail,  about  the 
motion  of  fluid  particles  through  interstitial  space.  These 
fluid  particles  are  identified  as  microscopic  differential 
volume  elements  each  assenblages  of  a  statistical  nuaber  of 
momentarily  contiguous  fluid  molecules)  to  tSiich  locations 
must  be  assigned  in  space  as  a  fuxiction  of  time  and  idiioh 
are  seen  to  trace  the  streamline  paths  between  sources  and 
sinks. 

We  find  it  instructive  to  describe  in  qualitatively 
the  tremendous  conqplexity  of  idiat  actually  happens  during 
flow  transfer  in  porous  media  systems,  for  out  of  tl:ls  we 
are  led  to  a  new  derivation  of  Darcy' s  law  from  Navier^tolces 
theory  which  avoids  the  suspect  features  of  previously 
offered  derivations  (Bidabert,  Philip, Ball,  Miller  and  Miller, 
etc.).  More  important,  our  an'-'.^ysis  has  the  advantage  of 
identifying,  a  priori,  systems  a  d  conditions  where  use  of 
Darcy' s  law  can  be  expected  to  fail  in  describing  the  xlow 

*This  manuscript  is  represented  as  being  a  collection  of  some 
ideas,  tentatively  taken  in  the  development  of  a  hew  approach 
to  porous  media  flow  problems.  A  large  amount  of  revision  and 
extension  is  anticipated  before  these  ideas  are  to  be  offered 
for  publication. 
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procMS  under  Invastlgation. 

Xn  the  end«  ee  wleh  to  treat  porous  aiodia  flow  systens 
according  to  the  aethods  of  potential  theory,  but  this 
representation  is  possible  only  to  the  extent  that  we  can 
assign  inacroscopically  neaningful  values  of  a  perneability 
tensor  to  volune  eleswnts  of  the  system  large  enough  to  c<»k- 
tain  a  "statistical”  asoount  of  pore  space,  but  small  with 
reference  to  the  relevant  gross  dimensions  of  the  systems  of 
Interest.  We  must  also  be  able  to  say  that  a  field  of  flux 
and  a  field  of  force  can  be  constructed,  %dklch  in  some  way 
approximate  a  one  to  one  correspondence  to  the  values  of 
velocity  vector  and  potential  of  the  microscopically  defined 
fluid  particles,  we  thus  construct  an  analog  apace,  each 
point  of  which  represents  discrete  volume  elements  in  real 
space,  and  we  regard  this  analog  space  as  a  continuous  media. 

We  find  that  a  continuous  media  representation  of  porous 
media  flows  cannot  be  undertaken  in  all  cases,  even  idien 
attention  is  limited  to  the  laminar  motion  of  incoaqpresslble 
fluids. 

Xn  the  microscopic  frame  of  reference,  there  will  be 
continuity  of  streamline  paths  whiure  these  will  be  defined 
as  the  trace  of  fluid  particles  (i.e.  a  microscopic  differential 
volume  element  which  is  an  assesblage  of  a  statistical  nusber 
of  momentarily  contiguous  fluid  molecules)  joining  sources 
to  sinks.  Xhe  paths  stay  not  be  fixed  in  space  with  time, 
but  intersections  leading  to  the  stagnation  points  characteristic 
of  certain  inviscid  fluid  flow  systems  will  not  be  encountered. 
The  fluid  particles  themselves  will,  in  general,  undergo 
successive  (and  perhaps  time«varient)  acceleration*^eceleratioas, 
leading  sometimes  to  a  dispersal  which  must  be  taken  into  account 
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if  BOlecular  diffusion  it  to  bt  repretontod.  To  tach 
fluid  particle  may  be  ascribed  a  certain  value  of  potential« 
Xinetic  and  pressure  energy,  as  «mi11  as  a  certain  value  of 
internal  energy  (including  surface  energy  if  the  fluid  particle 
is  located  at  an  interface).  In  the  limit,  the  pore  «fall 
boundaries  aure  the  loci  of  the  sero-velocity  streamline  paths, 
td)lch  is  to  say  that  the  adsorbed  fluid  particles  there  have 
sero  bine tic  energy. 

Continuous  matter  theory  deals  with  the  representation 
of  the  properties  of  matter  as  these  are  observed  in,  and 
defined  by  reference  to  asseaiblages  of  a  large  (statistically 
sufficient)  number  of  contiguous  molecules.  Such  assemblages 
themselves  can  be  quite  small,  as  illustrated  by  the  fact  that 
the  density  (or  other  property)  of  a  particular  material, 
held  every%ihere  at  the  same  temperature  and  pressure,  will 
be  the  same  as  a  sample  is  subdivided  (say)  below  the  limits 
of  visual  resolution.  The  subdivision  cannot  be  carried  with¬ 
out  limit,  however,  because  of  the  particulate  nature  of 
matter,  and  as  tho  volxame  element  dimensions  approach  molecular 
sires,  macroscopic  properties  of  matter  become  difficult  to 
define,  much  less  measure. 

In  this  essay  we  do  not  concern  ourselves  with  wbat 
gives  rise  to  particular  properties  of  matter,  that  is,  we 
make  no  appeal  to  kinetic  theories  nor  to  statistical  mechanics. 
We  choose  simply  to  postulate  that  we  have  a  continuous  medium, 
bounded  or  unbounded,  where  at  every  internal  point  a  property 
such  as  "demittance”*  can  be  defined  and  measured. 


*Demittance  refers  to  any  definable  and  measureable  property 
of  matter  whose  definition  and  measurement  are  left  unspecified 
in  a  particular  discussion,  to  avoid  the  introduction  of  any 
bias  which  would  come  from  the  discusser's  preconceiv^  notion 
about  the  property.  The  concept  and  the  term  were  prqpsed 
by  L.  deWitte,  Bugl^  Aircraft  Corporation,  1961. 


NETWORKS  REPRESENTING  POROUS  MEDIA  FLOW  SYSTEMS 

I.  FOUNDATIONS 

By 

Walter  Rose 


Introductory  Concepts 

In  the  end  we  shall  be  talking  exclusively  about  fluid 
flow  in  porous  media  systems g  but  much  of  what  will  be  said  has 
a  bearing  on  mass  and  energy  transport  phenomena  as  well  as  on 
momentum  transport;  therefore g  we  introduce  the  elementary 
concepts  in  a  quite  general  fashion. 

We  first  discuss  some  sort  of  a  continuous  conductor p 
without  specifying  whether  it  be  solid g  liquid  or  gaseous#  and 
whether  or  not  it  be  homogeneous  or  isotropic,,  If  conduction  is 
to  occur#  a  driving  force  must  be  acting#  whic.  is  to  say  that  a 
finite  energy  gradient  (e.g,  in  the-rmal  and/or  potential  and/or 
chemical  energy#  etc.)  will  be  found  directed  along  paths  joining 
the  sources  to  the  sink  boundaries;  moreover g  there  must  be  no 
zero“Conductivity  barrier (s)  completely  isolating  all  sources 
from  ail  sinks,.  The  conductor  may  be  thought  of  as  an  infinite 
medium  if  the  phenomena  of  interest  occur  (during  the  time 
interval  of  interest)  for  the  most  part  in  Internal  regions  far 
removed  from  the  influence  of  external  boundaries;  otherwise*  the 
external  geometry  must  be  a  specifiable  boundary  condition. 


Our  problem  is  to  seek  ways  to  describe  and  discuss 
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the  transport  phenomena,  given:  a  particular  conducting  medium 
having  a  particular  geometry,  a  particular  connection  to  sources 
and  sinks  of  given  strengths  an</^i  particular  initial  condition 
of  energy  distribution.  In  general  we  shall  find  it  necessary  to 
consider  unsteady  as  well  as  steady  states,  and  we  shall  have  to 
take  into  explicit  account  the  influences  of  the  inhomogenieties 
(which  say  that  contiguous  points  may  have  different  conductivities) 
and  the  anisotropies  (which  say  that  the  conductivity  at 
particular  points  may  depend  upon  direction^. 

Because  of  the  particulate  nature  of  matter,  it  is, 
in  the  final  analysis^  a  fiction  to  speak  of  a  continuous 
conductor  medium.  Oirb  the  other  hand,  as  long  as  one  is  dealing 
with  the  representation  of  the  properties  of  ma-ter  as  they  are 
observed  in,  and  defined  by  reference  to  assemblages  of  a  large 
(statistically  sufficient)  number  of  contiguous  molecules,  all 
functions  of  these  properties  may  be  thoi^ght  of  as  continuous 
throughout  the  mediumo  Therefore,  we  avoid  molecular  considerations 
by  taking  the  macroscopic  view  that  the  microscopic  volume  elements, 
each  comprised  of  a  large  number  of  molecules,  may  each  be  thought 
of  as  a  point  in  space.  The  infinite  ensemble  of  these  points 
(occupying  either  a  finite  or  an  infinite  space),  then,  is  our 
continuous  medium. 

Our  first  step,  therefore,  must  be  to  assign  values  of 
local  conductivity  to  each  point  in  space,  for  this  will  lead  to 
an  evaluation  of  the  effective  conductivity  of  the  medium,  and 
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to  a  prediction  of  what  transport  will  occur  for  given  initial 
and  boundary  conditions.  What  is  implied  by  this  statement  is 
that  we  have  a  partial  differential  equation  to  solve,  where  the 
dependent  variable(8)  of  interest  (flux,  potential)  are  functions 
of  both  time  and  position,  and  where  the  conductivity  of  the 
medium  (at  each  time  and  position)  explicitly  appears  as  a 
parameter o 

Indeed,  the  differential  equations  which  are  to  be 
solved,  in  the  most  general  cases,  are  nonlinear,  of  high  order, 
and  perhaps  inhomogeneous;  therefore,  great  problems  of  analysis 
are  met  with  from  the  start.  And  still  to  be  taken  into  considera¬ 
tion  are  the  added  complications  provided  by  the  stipulation  of 
unsymmetric  and/or  time-variant  boundary  conditions,  the  choice  of 
complex  initial  conditions,  and  the  circumstance  that  the  medium 
may  be  neither  homogeneous  nor  isotropic,'*’ 

We  choose  now  to  simplify  the  problem  by  saying  that 
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the  given  continuous  medium  can  be  subdivided  into  a  three- 
dimensional  array  of  a  finite  number  of  discrete  blocks.  We  then 
seek,  for  example  by  relaxation  methods,  to  describe  the  inter¬ 
action  between  each  contiguous  element.  In  effect,  we  hope  to 
come  to  numerical  results  by  the  traditional  procedure  of 
substituting  finite  differences  for  infinitesimal  differentials. 

The  suggestion  that  a  differential  equation  has  to  be  solved 
rests  on  the  assumption  that  a  satisfactory  basis  for  abstraction 
has  been  found  whereby  a  mathematical  model  (likely  more  imperfect 
than  perfect)  has  been  devised  to  represent  the  physical  system 
of  interest. 


and  by  discarding  the  notion  of  a  truly  continuous  conductor 
medium  in  favor  of  a  more  manageable  network  form.  That  is« 
each  block  of  the  network  will  have  lumped  in  it  the  appropriate 
characteristics  of  the  volume  element  it  represents  (  in 
accordance  with  some  averaging  process  still  to  be  defined) „  in 
a  manner  which  is  equivalent  to  saying  that  the  surfaces  which 
are  common  boundaries  to  contiguous  volume  elements  in  real 
space  are  represented  in  the  network  as  points  of  contact.  As 
a  consequence B  we  shal  1  reduce  the  degrees  of  freedom  of  the 
problem,  so  that  by  successive  approximation  we  hone  to  converge 
on  a  satisfactory  solution. 

To  restate  In  another  way  what  is  to  be  undertaken,  we 
are  saying  that  a  mathematical  point  in  real  space  is  a 
"physically"  infinitely  small  volume,  or,  as  described  by 
Landau  and  Llfshltz  (Fluid  Mechanics,  Volume  6  of  Course  in 
Theoretical  Physics,  Pergamon  Press,  1959,  Chapter  1) ,  a 
volume  element  particle  “...very  small  compared  with  the  body 
under  consideration,  but  large  compared  with  the  distances 
between  molecules...",  ^us,  we  may  expect  in  some  cases  that 
the  microscopic  differential  volvune  elements  roust  be  chosen 
say,  smaller  than  a  cubic  micron  in  order  to  think  of  a 
particular  medium  (say,  one  of  quite  small  dimensions  and  quite 
inhomogeneous)  as  a  continuous  conductor.  Alternately,  another 
(say,  larger  and  more  homogeneous)  particular  media  may  dis¬ 
play  all  the  characteristics  of  a  continous  conductor,  even 
when  volume  elements  having  dimensions  measured  by  orders  of 
cubic  meters  are  taken  as  the  "points"  in  mathematical  spaca. 

The  network  model  concept,  therefore,  is  derived  from  the  idea 
that  the  interaction  between  two  adjacent  collection  of  contiguous 

points  can  be  approximated  by  noting  the  Interaction  between 
two  other  jxsirts,  each  of  which  has  tlte  "average"  properties 
of  one  or  the  other  of  the  cc. .lections. 


The  Network  Model 


Let  us  first  address  ourselves  to  questions  of  potential 
flows.  During  the  laminar  transport  of  phases  (as  in  fluid  flov^„ 
components  of  phases  (as  in  diffusion)  p  o:c  energy  (e^g.  heat  or 
electricity)  through  conducting  media,  a  so-called  field  of  flow 
will  exist  which  will  be  representable  as  paths  joining  sources 
and  sinks  along  which  conduction  occurSo  The  trace  of  these 
paths  will  be  determined  by  the  inhomogeniety  of  the  medium^  by 
the  location  of  sources  and  sinks  and  by  the  geometry  of  the  other 
boundaries,  and  by  the  field  of  force  which  depicts  the  action 
of  the  force-intensity  vector..  If  anisotropies  exist  the  field  of 
flow  does  not  superimpose  on  the  field  of  force;  moreover!,  the 
flow  lines  may  deoend  on  the  initial  conditions  of  energy 
distribution  if  unsteady-stc 1 ensue,  .and  may  gradually  shift 
in  positiono 

In  order  to  anoroxim.ate  the  action  of  a  continuous 
conductor  by  a  network  modal,,  therefore,  one  n.ust  choose  a  form 
and  e  grid  spacing  such  th-at  the  flow  lines  of  the  network  are 
a  close  ('i:ep-wisc)  r apresent -.tion  of  the  flow  lines  of  the 
prototype.  Tnerefore,  as  a  rule  of  chv.iri’o,,  coarse  spaclngs  of 
network  elements  will  a  reasonable  aoproxiination  when  the 
medium  being  re'pres.;ntcd  is  ’  omoqeneous  and  i.sotropiCt,  and  for 
when  attention  is  limited  to  steady-states,  and/or  when  the 
boundaries  are  su.’h  iha  L  o  .-xr  r.'e  netv.nr'.c  form  can  be  constructed 
which  has  the  cenf igur-tlo-i  c.'oscj.y  nupt'rimposing  on  the  actual 


field  of  flow;  otherwise  fine  spacings  are  required„  which  is 
to  say  that  the  model  rust  represent  a  high  degree  of  subdivision 
of  the  continuous  conductor  media  prototype. 

The  foregoing  statementSf  of  course,  suffer  in 
exactness  from  their  own  relativityp  but  the  point  is  clear  that 
in  many  oractical  cases  it  will  be  possible  to  assign  upper  and 
lower  limits  to  the  number  of  network  elements,  optimized  so  that 
resolution  and  convergence  on  the  sought  answer  will  be  maximized, 
while  the  computational  labor  will  be  minimized.  Indeed,  it  is 
the  major  purpose  of  our  discussion  to  examine  how  to  arrive  at 
a  decision  about  network  size  and  form,  and  to  consider  how  much 
will  be  lost  in  accuracy  of  representation  in  the  effort  to  avoid 
an  unreasonable  anount  of  numerical  analysiso 

Figure  i  shows  one  plane  of  a  cubic  network  form."^ 

Each  contiguous  node  (grid  point,  mesh  point,  etc,)  is  to  be 
thought  of  as  representing  the  center  of  volume  of  contiguous 
macroscopic  volume  elements  (A^)  prototype;  therefore, 

the  number  of  nodes,  N,  is  a  measure  of  the  degree  of  subdivision 
which  is  intended.  The  nature  and  the  size  of  the  elements 
surrounding  each  node,  then,  are  taken  to  represent  the  relevant 
average  properties  of  the  protot/pe  medium  in  the  volume  element 
which  each  node  point  depicts. 


A  wide  variety  of  two  and  three  dimensional  network  forms  may 
be  imagined,  which  can  be  regular  or  irregular.  In  particular 
cases,  however,  the  geometries  may  be  such  that  it  mav  prove 
sensible  to  choose  a  particular  network  form  to  maximize  the 
overlap  of  the  flow  paths  in  the  model  and  the  prototype  systems 
for  a  given  value  of  N.  On  the  other  hand,  we  shall  limit 
reference  to  cubic  forms  as  they  are  easy  to  discuss,  and  in  fact 
always  can  be  used  if  N  is  large,. 


In  effects  we  started  with  a  notion  of  a  continuous 
inediunip  and  took  the  (usually)  safe  position  that  a  statistical 
collection  of  molecules  is  "small”  enough  to  be  representable 
as  a  mathematical  point.  The  sense  of  constructing  a  network 
model  is  that  v/e  now  have  taken  a  statistical  collection  of 
contiguous  points  (i.e.  the  ricroscopic  differential  volume 
elements;.  dV,  each  of  which  contains  a  statistical  number  of 
molecules) e  and  we  call  this  macroscopic  differential  volume 
elements  also  a  point  (i.e.  a  node  point)  in  the  network 

model  space.  Thus,  we  have  constructed,  in  network  form,  a  new 
quasi-continuous  medium,  the  contiguous  points  of  which  represent 
a  sensibly  finite  distance  in  real  space. 

It  follows,  that  in  transforming  from  real  soace  to 
network  mode]  space,,  we  gain  little  or  nothing  in  the  analytic 
formulation  of  the  problem,  even  though  the  exact  initial  and 
boi.'ndary  conditions  and  the  exact  local  conductivity  of  the 
prototype  system,  in  general,  are  simplified  (i.e.  roughly 
aoproximated)  in  the  network  model  analog."'’  Thus,  Figure  2 
depicts  how  a  given  network  form  and  spacing  can  only  anproximate, 
but  not  match  the  conditions  and  performance  of  a  particular 
conduction  system,'*"*'  But  in  soite  of  the  loss  of  resolution  when 
the  problem  is  cast  into  its  network  analog  representation,  the 
exact  solution  of  the  network  model  problem  involves  a  similar 

+  Elsewhere  we  discuss  at  length  how,  in  the  case  of  fluid  flow 
through  porous  media,  a  simplicity  is  introduced  in  the  fact  that 
a  more  tractable  law  of  force  applies  to  the  network  model  system 
^.e.  Darcy’s  law)  than  to  the  microscopic  prototype  system 
{l,e.  the  Navier-Stokes  set  of  equations), 

++  The  prototype  flow  path  of  Figure  2  meanders  in  a  way  related  to 
the  inhomogeniety  and  anisctrooy  of  the  svstem,  as  well  as  to  the 

irr:posed  and  'uunciiir'-  -d;  c.i oi  s. 
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analytic  vormulation  as  characterizes  the  solution  of  the 
prototype  problem..  In  both  cases,  conservation  of  mass  and 
energy  must  be  established  for  every  point  in  the  system  (i.e.  the 
continuity  statements),  and  an  equivalent  law  of  motion  must  be 
stated  and  applied  to  every  point  in  the  system.  This  means,  in 
addition,  that  the  conductivity  tensor  must  be  known  at  every 
point,  whether  we  are  talking  about  the  prototype  system  or  the 
network  analog  thereof. 

On  the  other  hand,  and  as  to  be  expected,  advantages 
do  result  by  transforming  a  conduction  oroblem  into  its  network 
representation.  In  the  first  place,  if  the  prototype  system  is 
characterized  by  a  high  degree  of  inhomoqeniety,  explicit  solutions 
are  not  to  be  expected  even  for  cases  of  simple  initial  and 
boundary  conditions.  Similarly,  if  the  prototype  is  quite 
homogeneous^  (whether  locally  there  is  isotropy  or  not), 
explicit  solutions  are  available,  it  at  all,  only  for  cases  of 
simple  initial  and  boundary  conditions.  It  will  be  seen,  therefore, 
that  the  gain  achievect  by  reducing  the  problem  to  its  network 
model  form  is  three-fold.  If  an  analytic  solution  is  possible  at 
all,  one  immediately  gains  the  advantage  that  conditions  and 


We  define  homogenietv  to  describe  a  medium  where  the  conductivity 
tensor  at  all  contiguous  points  will  be  the  same,  or  nearly  so, 
and  will  have  the  same  oonf In  consequence,  many  types  of 
inhomogeniety  (i.e,  heterogeniety)  may  be  imagined,  including 
bedded  media,  graded  media,  media  where  the  inhomogeniety  comes 
only  from  the  fact  that  there  is  an  irregular  orientation  of 
anisotropic  contiguous  elements,  etc.  In  any  case,  it  is  clear 
that  any  network  model  will  tend  to  be  less  inhomogeneous  than 
the  prototype  from  which  it  was  derived,  and  will  tend  to  have 
an  altered  anisotropy  character.  In  a  corresponding  sense,  it 
is  to  be  noted  that,  in  the  limit,  a  medium  which  is  completely 
inhomogeneous  in  the  fri cros ;pnic  pensp  (by  the  de'f’n.nition  given 
abov?)  would  have  th^  c.-two.'-d  'r:a  ,,  p  ^ 
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characteristics  need  be  specified  only  at  a  limited  number  of 
points  (ioSo  at  the  N  nodes),  instead  of  continuously  through¬ 
out  space.  Secondly,  if  the  number  of  nodes  is  not  too  large, 
closed  (implicit)  solutions,  involving  the  simultaneous  solution 
of  some  M  equations  with  N  unknowns,  are  practical  to 
employ.  Lastly,,  network  problems  are  by  their  formulation 
already  reduced  to  the  point  where  standard  relaxation  and  other 
(explicit)  numerical  methods  (cf.  for  example,  the  Hypercircle 
method  of  Synge)  cnn  be  applied. 

Network  Model  of  Continuous  Media 

We  choose  now  to  Illustrate  some  points  by  making  refer¬ 
ence  to  a  three  dimensional  heat  flow  system,  for  example,  that 
having  the  geometry  of  Figure  3,  and  made  up  of  a  large 
(essentially  Infinite)  number  of  single  crystals  of  various 
atomic  species.  Each  crystal  will  be  thought  of  as  being  at 
least  as  large  as  whatever  is  required  so  that  its  lattice  con¬ 
tains  the  (so-called)  statistical  number  of  atoms;  therefore, 
if  this  mixture  of  crystals  is  pressed  together  to  form  a  solid 
mass,  we  can  say  we  are  dealing  with  a  (virtual)  continuous 
medium.  Further,  we  shall  say  that  the  medium  is  inhomogeneous, 
at  least  on  the  microscopic  scale,  because  we  postulate  that 
each  species  of  crystal  has  its  own  value  of  specific  heat 
conductivity,  (Tij  .  To  be  general,  we  say  that  some  of  the 
crystals  undergoing  radio-active  decay  (say,  of  sensibly  long 
half-life)  so  that  we  aay  of  then  as  internally  distriboted 

haat  sources.  Other  (aocteznal)  sources  and  sinks  of  given 
strengths  are  shown  in  Figure  3A,  together  with  tiie  trace  of 
certain  arbitrary  flow  linas  of  intereat. 
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We  do  not  let  transfer  of  energy  by  convection  or  radiation 
enter  into  the  formulation  of  this  problem.  1herefore« 
to  describe  heat  flow  (by  conduction)  under  these  conditions* 
we  need  to  know  the  Initial  distribution  of  temperature  at 
every  (crystal)  point  in  the  system*  and  the  strength  (and 
its  time  variation*  if  any)  of  every  source  and  sink;  moreover* 
we  must  know  the  thermal  conductivity  tensor*  of  each  single 
crystal  point*  and  %fe  must  have  the  geometry  of  all  surfaces 
specified  across  which  there  is  no  flow  (i.e.  the  luqpermeable 
boundaries) . 

Now*  for  each  point  in  the  system,  except  on  the  source 
and  sink  boundaries  and  the  inteamally  distributed  source 
points  *  a  continuity  statement  can  be  applied  to  the  effect 
that  the  divergence  of  heat  flux  is  proportional  to  the  Change 
in  heat  content  with  time*  and  Fourier*  s  law  of  heat  conduction 
can  be  applied  to  give  as  the  govemlng  differential  eguation 
of  motion: 

where t 


•  Mioroeooplo  flux  per  \mit  curea  vector, 

■  microscopic  value  of  potential  at  (x,y,a,t), 
^  microaoopio  value  of  local  capacitance, 

«  local  value  of  conduotlvlty  tenaor. 


TO  be  noted  is  the  fact  that  Equation  1  is  most  readily 
integrkble  when  steady-at^tes  aura  being  discussed  (i.e.  the 
right-hauAd  meniber  is  e<ji.r  >.  to  zero*  so  that  Equation  1  reduces 


*It  is  a  consequence  of  continuous  media  theory  that  a  symmetric* 
second-rank  tensor  will  desc.^ilM  the  most  general  case  of  aniso¬ 
tropy  (cf.  transmission  of  heat*  electricity*  electromagnetic 
radiation*  acoustic  energy*  fluid  flow*  diffusion*  etc.  through 
amisotropic  media) . 


to  the  Laplaclan  fona) ;  but  even  with  this  ■inqplification,  the 
io^sltion  of  aaynetrlc  axul/or  non-linear  initial  and  boundary 
condition  in  moat  caaea  will  introduce  great  analytic  difficultiea 

One  device  for  circumventing  analyala  probleaia,  and  in 
fact  the  one  «diich  we  propoae,  le  to  conatruct  an  analog  net¬ 
work  of  the  given  continuoua  conductor  ayaten,  aa  ahown  in 
Figures  3Bj  3C  xmd  3D.  (In  ways  to  be  deacribed  later,  we 
can  determine  that  Figure  3C  repreaenta  the  optimum  model  of 
the  prototype  system,  as  the  coarse  grid  sought,  and  the 
fine  grid  of  Figure  30  would  involve  too  much  coaptation) . 

These  models  ho%rever,  are  at  best  only  approximately  equivalent 
to  Figure  3 A. 

In  any  case,  the  network  model  of  Figure  3C,  in  turn, 
may  or  auiy  not  be  analyzable  by  exact  methods,  largely 
depending  on  the  ntudxer  of  node  points,  N,  involved,  and  to 
a  leaser  extent  depending  on  symnetry  and  homogeniety  con¬ 
siderations.  In  order  to  anticipate  the  difficultiea  of  the 
most  general  case,  however,  let  us  assume  that  analytic 
solutions  are  impractical*  or  iaqposaible. 

we  now  assume  that  a  flow  equation  analogous  to  Equation 
1  can  be  stated  which  will  apply  to  the  node  points  of  the 
network  model.  Since  the  latter  arbitrarily  has  been  put  in 
a  cubic  form,  %fe  choose  the  coordinate  axes  to  correspond  to 
the  principal  directions  of  the  network,  giviz^: 


0 


*Our  high  speed,  digital  ,conpter,  ILLIAC,  in  general  cannot 
handle  more  than  140  simultanecjus  equations  of  140  unknowns 
This  refers  to  linear  polynomial  matrices. 


(12 


It  Iri  thua  Inpliad  that  the  terms  of  Equation  2  are  to 
be  thought  of  as  averages  obtained  from  the  terms  of 
Ecuation  1  according  to  the  definitions: 


T(p'  -  I'-  r 


J%i/ 


#  -  J -14 

^  “  v{"  ^  ^ 

Now  it  will  be  clear  that,  ais  N  approaches  infinity. 


Equations  2  will  give  the  same  result  as  would  be  obtained 
by  an  analytic  solution  (inhether  kno'im  or  not)  of  Equations  1, 
at  least  to  the  extent  that  the  successive  approximations 
taken  in  connection  with  So:uation  2  leade  to  sensible  convergence. 
Therefore,  the  question  at  hand  is  t/nat  value  of  N  shall  be 
chosen  which  will  give  the  numerical  accuracy  6f  results 
which  is  wanted  without  entailing  an  vmreasonSble  amoxant  of 
computational  labor? 

Several  reasons,  arrived  at  intuitively,  can  be  cited 
why  the  network  form  o£  Figure  3B  is  less  satisfactory  than 
that  of  ?igvire  3C  in  representing  the  prototype  system  under 
discussion  (cf .  Figure  3h) .  We  wish  to  achieve  a  high  degree 
of  superposition  between  the  fields  of  force  and  flux  in  the 
phototype  and  model  sy, it  tens. 


Fr.GURES  1  and  3  ’.vi:,  1  be  supplju^d  a"!*,  a  latar  date 


1  Iroto  ’r,;/.  --  real  (micro&ooplc)  ypaco,  v/hore  wixhin  Vie  reciona 

A  and  B,  fluid  pf.ad.j,cies  have  local  end  'c-iae-dependunt  values 
of  poie;ixi;.il,  f.* U'-  cl  of  th'?  ve?.ocity  vector,  v* .  , 

Li  I  1  oa';  !  “  i^piasi-GOUtinuouE  space  uaae  up  of  contiguous  macro- 
3>- opic  luaps  or  real  space 

HI  Network  model  form  of  II » 


yiGiriii 
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GBMBRALZZBD  fLOUT  PKOBUM  XH  POROUS  lODXA 
lay  Bao  Chanaapragada 


Tha  tMo-dlaMnalonal  flow  pxOblcai,  vith  a  danuMurSbla 

nuaiber  of  arbitrary  abapad  gralna  in  a  pora  cbannal,  ia  VMlysad 

analogous  to  the  taraatnent  of  P.  B.  Richards  and  sesMi  applica<- 

(2\ 

tions  of  tha  author*  s  doctoral  thesis '  have  been  introduced. 

Aa  pointed  out  by  Richards  each  of  tha  contours  of  tha 

grains  can  be  replaced  by  a  sheet  of  vortasc-like  singularities 

such  that  tha  velocity  induced  by  each  aleawnt  of  tha  sheet 

satisfies  the  bouzidary  conditions  of  tha  flow  on  tha  channel  %falls. 

Tha  cosqplax  velocity  on  the  contours  satisfies  tha  singular 

integral  aquation  with  Cauchy-typa  kernel  taken  over  tha  contours 

alone.  This  can  be  reduced  to  a  linear  Fredholm  integral 

equation  of  the  second  kind  «^th  a  continued  kernel  for  the 

tangential  velocity  on  the  contour  and  the  solution  of  the 

(3) 

Rredholm  equations  follows  naturally  •  The  method  of  attack  of 
the  flow  problem  using  singular  integral  equations  offers  definite 
computational  advantages  and  it  suffices  to  quote  Richards « 
■Previous  analytical  work  on  channel  flow  has  been  more  restricted 
in  its  applications  and/or  complicated  computationally.* 

With  no  loss  of  generality  let  the  pore  channel  be  an 
infinite  strip  of  width  r  and  free 

stream  velocity  of  the  flow  in  the  positive  x<direction  be 
V  •  (v«l«2>...n)  is  the  complex  velocity 

by  the  entire  vortex  sheet  coinciding  with  the  oontours 
of  the  grains  C^. 


iMt  £(s)  1M  th»  coapl«9c  pottttial  of  tbo  iroxtOK-Ilko 
■ingularlty  Xocatod  at  oa  C^« 

Wa  now  Betlk  an  analytic  aoIutlMi  to  the  ehaanol  flow 
In  the  fon 


«(s«r)  «  V  +  U^Ca.r) 


V  (acM  +^)  m  0 
V  * 


(1) 


The  Induced  velocity  corresponding  to  C 


”»“>'>  -  air 


^  h(t. 


r)  ^  (z«t,r)  dt 


(2) 


where  the  cooplex  potential  of  a  unit  vortex  at  t  d  la  given 
by 


(•n'yr) 

esq^  (TTz/r)  +  esqp  (TTC/r) 


(3) 


and  the  coaplex  velocity 


SOL^  i  ^<t,z.r) 
ds  ^  t-e 


(4) 


where 

u(t,s,r) 

Therefore 


U^(z,r) 


T  r  t  >  a  i.r,,.! _ _ 

"  r  I  eaqp ( Tf/r )  (t'^z )  -  1  exp(TT7z)  (1>«)  ♦  1 
1  AY  h(t,r)  u  (t,s,r) 

““ilrr  T  - 


(5) 


t  -  z 


(6) 


In  aqiuition  (1)  thm  ooiitdltion  m  o  on 


iapllMi 


haa  juap  valuoo  at  C  ainoa  V  la  a  ccnatant. 

V 


in 


tJ^(a,r) 


by  attlelny 


Wa  solva  for  thaaa  Juaqp  valuai 
uaa  of  the  wall-loaoiai  Cauchy  Intagrala* 

Por  the  PiaaMlJ  foxanXaa  to  hold  alony  c^$  wa  have  the 
following  hypothaaiai  Zf  the  danaity  function  g(t«8)  aatiaflaa 
a  Lipaehita  c<»iditias  in  a#  in  the  neighborhood  of  a  anooth 

and  alao  aatiaflaa  a  Lipaehita  condition  la 

of  the  Cauc)^  integral 


Jordan  curve 
t  on 


C^,  then  the  principal  value 


G(a) 


aiTi 


f 


o(to) 


g(t.«) 

t-a 


dt 


(7) 


exiata  for  ^  C  ^  well  aa  the  limiting  valuea  6  (t^) 

and  6*(t  )  aa  the  point  t  la  approached  from  the  left  (’!>) 

of  and  the  right  (•>)  of  C^. 

Ttaua  the  Plemelj  fOnmlas  can  be  nrritten  aa 


G^(t^)  -  o"(t^)  - 


(8) 


Applying  theae  relations  to  the  induced  velocity  U^(8,r)  which 
belongs  to,  the  class  of  Cauchy  integrals,  we  Obtain  an  integral 
equation  of  the  second  kind. 

Assuming  that  the  vortex  strength  k(t,r)  satisfies  a 
Lipaehita  condition  on  C^,  we  have  from  the  Plemelj  relations 


^  “  20(t^.r) 

U^(t^,r)  -  H  +  n/(t^,r)  ^ 


(9) 


(10) 


I 


I 

I 

I 


-4.- 


and  £xo«  (1)  «•  h«v« 

T>^(t^,r)  ■  V  + 

-  V  +  U*(t  ,r) 
o  o 


But  vi^it.i')  »  0  Inside  C 

O  V 


From 


- 


(i),  (13), 
U(t^,r)  « 


.  -v 

2 

(9)  and  (11)  we  have 

2 


lc(t  ,r)  tt  (t  .t  .r)  =TJ^(t  ,r) 
o  o  o  o 


(11) 


(U) 


(13) 


(14) 


Alt  u(t  ,t  .r)  «  1 
Q  O 

Hence  the  vortex  strength  is  equal  to  the  tangential  velocity 
along  C^,  i.e. 

k(to2r)  »  lF(t^,r) 


■5' 


Tims  upon  sulMitltatioB  in 


ll^(s#r)  ■  ?V 


(6) 


«s  olstslB  tlis  latsgrsl  sqpistlsn 

Vf-ttidJ  ait 

t  -  s 


CIS) 


%flmrs  tbs  latsgrsl  oa  tbs  sigbt-bsad  slds  of  sgostlOB  (IS)  is  tbs 
prlBoips7>  valus.  Thus  if  tbs  chaansl  flow  paroblsa  with  gssia 
bouadariss  adaits  a  solution  socb  that  vsloeity  ll^(t«r)  at 
satisflss  a  Lipschits  conditioa  tbsn  it  is  a  solution  of  tbs 


intsgral  squatlon  (15)  and  tbs  vsloeity  fisld  is  rsprsssntsd  by 


0*(s«r)  «  V 


ir(t.r)  u  (t.s.r) 
t  -  z 


(16) 


If  the  do  iatsrsectf  l.e.  i£  tbs  grain  bouadariss  ars  ia 
contact*  than  a  study  o£  tbs  bShavior  of  tbs  singular  intsgral 
sqpiation  (16)  at  tbs  points  of  intsrssction  o£  is  to  be  nads. 

Tbs  analysis  for  tbs  point  of  intsrssetion  t*  ( 4  Cy) 
is  analogous  to  tbs  above  discussion  with  tbs  following  nodiflsd 
PlasnlJ  fomulas  ; 


oJ(t*)  +  0^'^'*'^W)  -  2Q(t*) 


(t*)  -  G^^^^^(t*)  •  g(t*,t*) 


Riehaxds,  *R«prM6ntation  by  vort«x*liX«  slagnlaritiMi 
of  obataeXo  flow  in  m  straight  ohaaBsl'**  tXAM  Vbl*  2, 
Janm  1961 • 

Chaanapragada#  "Singular  iatsgral  aguatioos  with  an 
Infinity  of  intaxsaetlng  area**  9h.O.  tbaaia*  IXaivascaity 
of  Ulinoio*  1961. 

TTjitsinaky*  "Singular  integral  eguationa  with  Cauclqr 
kernels**  Trans.  American  Mathasmtlcal  Society*  >A>1.  60* 
No.  2*  Septesiber  1946. 

Miskhelishvlli*  "Singular  integral  aguatioos** 

P.  Noordhoff  N.V.*  Groningen*  Bollaad*  1953. 
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?L^»m.1:iv  8iippl«aeotary  »ta 

Final  Import 

By 

Rao  Chaimapragada 

Stumaary  of  Anoy  Chaalcal  Mfarfara  Projact  CNXr»517 
"Analytical  Study  of  the  Behavior  of  Fluids  in  Porous  Solid 
Media". 

author  was  officially  appointed  to  the  contract  as 
a  jresearch  associate  from  Febzruary  1>  1960  to  carry  out  the 
analytical  phase  of  the  study. 

mie  following  is  an  itemised  account  of  the  progress 
made  by  the  author  during  the  period  of  his  appointment. 

I} .  Some  basic  concepts  and  definitions  with  regard  to  the 
geometric  properties  of  minimal  surfaces  and  constant 
curvature  surfaces  were  defined. 

2)  .  In  the  study  of  immiscible  flow  through  porous  media  a 

general  definition  of  capillary  pressure  and  the  eguations 

(Ij 

of  motionCfron  microscopic  viewpoint)  were  arrived  at. 

3)  .  During  the  study  of  surfaces  of  discontinuity  the  writer 

was  led  to  the  study  of  singular  integral  equations. 

■This  resulted  in  the  cospletion  of  the  doctoral  thesis 
(Z) 

of  the  author  .  The  thesis  however  was  purely  from 
the  mathematical  standpoint  as  it  vas  submitted  in  the 
Department  of  Mathematics. 

4)  o  An  application  of  the  author' s  thesis  has  been  indicated 

(31) 

in  the  report  along  the  linea  of  P.  B.  Richards. 

5)  .  An  analysis  of  the  flow  problem  from  the  fluid  dynamic 

standpoint  was  initiated  and  a  qualitative  analysis  is 
included  .  T^ie  analysis  can  be  further  improved  from 


(2) 


th«  aathMMtieal  standpoint.  Bowsvar*  wich  tba 

aid  o£  tha  physical  pictura  as  prasantad  in  tha 
14) 

raport'  \  it  is  baliavad  that  tha  sttthasiatieal 
analysis  will  now  ba  facilitatad  to  a  graat  axtant. 

Racamandatlons  for  furthar  study: 

a)  Zt  is  strongly  advisabla  to  confin  axpari- 

nantally  the  actual  flow  pattern  of  the 

streamlines  and  the  prasanca  of  tha  moving 

(4) 

vortices  as  predicted  by  tha  analysis  . 
b}  There  is  a  singularity  at  a  point  where  tha 

interface  neats  tha  wall  of  the  capillary  tube. 
If  tha  moving  interface  does  not  change  its 
shape  tha  boundary  condition  at  tba  point  of 
contact  of  tha  interface  and  tha  wall  is  to  ba 
Justified  since  there  is  a  sudden  jusqp  in  the 
magnitude  of  tha  velocity  (from  saro  to  a  finite 
value) . 


(3) 


■  ^  -t  ^  J  ^  ' 


X).  "CapilXazy  Praasura  and  aurfaca  diacMtiattl^  In  Poxoua 
Madia*  by  Mao  duumapragada  and  Waltar  Moaa»  Jour, 
of  CkiqpliyaicaX  Maaaaroh«  April  1961.  Vol.  66  #4. 


2).  "Singular  intagral  aquationa  with  an  infinity  of 

intaraacting  araaa*  by  Mao  Cbannapragada.  Doctoral 
tbaaia,  Oapartaant  of  MathaBAtica*  Dnivarai^  of 
Zllinoia^  Juna  1961. 


3) .  *Oanaraliaad  flow  in  Poroua  Nadia”  By  Mao  Cbannapragada. 

(Saa  pravioua  aactlon) . 

4)  .  "TMo  diaAnaional  Analyaia  of  an  iaaBiaelbla  flow 

through  a  capillary  tidia*  by  Mao  Cbannapragada. 

CSaa  naxt  aactlon) . 
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TND  01MB»ZC»ME.  MIAUSrSZS  OT  Ml  XMIZ8CZBLI  WUm 
TBBOm  A  CAPZLlAKr  fOBI 

By 

Raci  Channapragada 


Lat  tha  eapillaxy  tuba  ba  plaead  boricontallp  and  tha 
prasauaca  gradiaat  ba  tha  torca  eaualag  tha  flow  of  two  laalaolbla 
flulda  acroaa  tha  tuba.  Tha  flow  la  aaauBMd  to  ba  ineon* 
praaalbla  vlacoua  and  laminar.  It  la  wall  known  that  tha 
Intarfaca  (aurfaea  of  diacontinulty  batwaan  tha  two  flulda) 
moving  acTOaa  tha  tuba  aatlaflaa  Laplacaa  condition 


%diara  R^ 


and  R2  ara  principal  radii  of  curvatura  and  K  la 


aaauBMd  to  ba  a  conatant.  During  tha  flow  It  la  aaaumad  that 


tha  Intarfaca  doaa  not  change  ahapa  onca  tha  flow  haa  boon 


atablllaad  (tha  praaaura  dlffaranca  acroaa  tha  tid>a  la  kept 


conatant) .  Ihla  haa  boon  verified  akparimantally  by  Roaa 


fOr  laminar  flow  tha  velocity  profile  (Figure  1)  la 
parabolic  and  at  tha  Intarfaca  Itaalf  tha  valocl^  at  each 
point  of  tha  Intarfaca  la  tha  aama. 

lha  quaatlon  that  naturally  arlaaa  la  aa  to  how  doaa 
tha  parabolic  velocity  pcroflla  tranaform  Itaalf  to  that  of 
tha  Intarfaca  (Figure  1). 

From  tha  pbotographa  of  Foattlngar^^^it  la  daduead  tttat 
tha  advaraa  praaaura  gradient  together  with  friction  near 
tha  wall  dataalna  tha  proeaaa  of  a^aratlon  (Figure  2)  • 
lha  praaanea  of  a  thin  plate  along  tha  plana  of  ayMMtry  at 
right  anglaa  to  tha  first  wall  (rafaranca  3,  Figure  2.14 
and  2.15)  gives  rise  to  tha  two  symmatrically  placed  vortaclas. 
Mow  Imagine  tha  same  vertical  wall  with  thin  plate  to  ba  set 


(2) 


la  notion  teeauso  of  tho  pcooonro  foceo  aotlag  on  tho  wall. 

Zt  can  again  bo  ahoioi  aaqporlaantally  that  tboro  vlll  atlXl  bo 
tMD  stablllsod  vortoelos  aovlag  with  tho  wall  (slallar  to 
Flgoro  2) .  Tho  atrmgtb  of  thooo  vortoeloo  vlll  bo  lowar  tflMa 
In  tho  caoo  of  tho  otatlonnxy  nail  boeanoo  of  tho  doeroaoo  In 
tho  xolatlvo  voloelty  of  tho  flow  Juot  ahoad  of  tho  vortical 
wall. 

Zf  Inotoad  of  tho  thin  plato  at  tho  cantor  (Figure  2) 
two  and  platos  are  placed  cm  tho  ends  of  tho  vortical  %nill> 
we  will  oboorvo  a  fountain  effect  (Figure  2a)  Ohoro  tho  flow 
rovoraos  at  tho  comers  and  tho  vortoelos  are  again  synsiotrlcal 
but  now  displaced  from  tho  cantor  of  the  novlng  wall  to  tho 
comers. 

A  slnllar  phonoaona  takes  place  In  tho  case  of  tho  amvlng 
intorfaco.  Boro  tho  vortical  wall  Is  replaced  by  tho  Inter¬ 
face  and  tho  cmd  platos  by  the  walls  of  tho  capillary  tid>o. 

Consider  tho  region  A-B  of  Figure  3.  The  stroan  linos 
aoparato  froai  tho  relative  stagnation  point  (novlng  with  tho 
Intorfaco)  and  tno  stable  vortoelos  are  fomad  novlng  with 
tho  Intorfaoo  (Figure  3) .  Also  a  secondary  flow  Is  sot  up 
duo  to  tho  fomatlon  of  tho  vortoelos  causing  tho  fluid 
particles  near  tho  wall  to  bo  pushed  baokwords. 

Iho  fluid  being  viscous  and  for  low  Boynolds  nusbors 
tho  secondary  flow  effect  cannot  bo  carried  too  far  fron  tho 
Interface  as  shown  schanatloally  by  tho  velocity  profiles  In 
Flgurf  4.  Iho  regions  of  secondary  and  prlnary  flows  are 
soparatod  by  tho  Uniting  linos  as  ahOMu  la  Figure  5. 

BOW  consider  tho  volocl^  profiles  across  tho  Intorfaoo. 

As  tho  Interface  Is  novlng  wltti  a  constant  velocity^  the 
volocl^  vectors  on  either  side  of  tho  Intorfaoo  are  equal 


(3) 


(in  aganitud*  and  diraotlea) ,  (valeoi^  profllaa  and  P^) 
(Figura  6) .  ibat  ia  tha  flow  on  aitbar  alda  in  tha  lanadlata 
viciai^  of  Idia  latarfaoa  ia  in  tha  aaaw  diraetion*  Xn  tha 
ragion  ac  of  Pigura  3  idiara  ia  an  inraraa  fountain  affaet  in 
tha  atnaaai  lina  flow.  Tbia  ia  again  oauaad  by  tha  praaanoa 
of  tha  two  vortaeiaa  aa  ahown  in  Pigura  3. 

Iba  atagnation  point  for  tha  ragion  BC  oecura  at  a  point 
hatwaan  tha  vortaeiaa  and  just  hafOra  tha  coanMoeanant  of  tha 
invaraa  fountain  flow«  indioatad  by  tha  daahad  linaa  aaparating 
tha  aacondary  flow  ragion  froai  tha  laain  flow. 

Sbia  analyaia  givaa  us  a  physical  pictura  of  tha  flow 
phaaoaiaaa  ia  tha  capillary  tuha.  lha  Sbova  analysis  has 
haan  earriad  out  strictly  fron  physical  reasoning  and  it  is 
hoped  that  a  aukthasMtical  analysis  of  tha  above  would  rain- 
force  tha  above  analysis. 

Tha  idea  of  tha  fountain  affect  was  originated  froai  tha 
discussions  trith  Professor  w.  O.  Rosa  and  tha  photographs  of 
Poattingar  gave  an  insight  into  tha  physical  reasoning  of  tha 
above  analysis. 
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ABSTRACT 

If  on«  end  of  a  capillary  filled  with  air  or  any  vapor 
ie  placed  in  contact  with  a  reoervpir  of  wetting  liquid* 
flow  into  the  tube  will  be  initiated  spontaneously  and  isi* 
pulsively.  An  approxisuite  analysis  of  the  subsequent  motion 
considers  the  inibibinyt  column  as  a  single  (increasing) 
mass  and  yields  an  oaniinairy  differential  equation  for  the 
extent  of  rise  as  a  function  of  time.  Ibe  setting  of  the 
Initial  conditions  (the  velocity  and  position  of  the  interface 
between  the  liquid  and  vapor)  presents  some  difficulties: 
earlier  treatments  produce  scxae  inconsistencies  and  the  re¬ 
sulting  integrals  do  not  exhibit  the  initial  impulsive 
transient.  In  this  paper  general  initial  conditions  «diich  are 
physically  and  mathematically  consistent  are  developed.  The 
duration  of  the  transient  is  given  as  a  function  of  the 
initial  conditions*  and  of  the  starface  tmnsion*  viscosity 
and  density  of  the  liquid.  An  upper  boaand  for  the  maximum 
velocity  is  derived. 
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ZHTMODUCTZOR 


An  incoiBpres8ibl«  wsttlng  liquid  of  donaity  p  and 

viscosity  1.  axid  driven  toy  capillary  forces  is  aivanetng 

2 

into  a  capillary  of  unifom  circular  cross-oection  )f  r 
and  length  and  displacing  air  or  aiqr  vapor  of  lov  density 
(  p  )  and  low  viscosity,  (  )  •  Bven  for  this  slaqple 

exaa^le  of  cbplllary  flow,  an  exact  oiathenatical  analysis 
%iK>uld  involve  the  solution  of  an  intractatole  two-region 
(liquid  and  vapor)  toovindary  value  problem  consisting  of  three 
partial  differential  equations  and  appropriate  bovu&dary 
ccmditions.  [[l3  Further  coaqplications  arise  because  the 
capillary  driving  force  is  a  function  of  the  norentary  fluid- 
fluid  interface  geometry,  and  thus  nay  be  time-dependent; 
and  even  if  conqpreasltoility  is  neglected,  it  is  difficult  to 
set  well-defined  conditions  at  this  surface  since  a  model  for 
this  third  '  region*  is  required  £,2,  $3  .  An  approximate 

treatment  considers  the  inbltoing  liquid  column  as  a  single 
(increasing)  mass  subject  to  driving  and  retarding  forces; 
the  viscosity  and  density  of  the  viqpor  are  usually  neglected. 

This  approach  goes  back  at  least  to  HOlf f  4  and 

has  been  adopted  toy  sumy  others  *  C)  *  yields  an  ordinary, 
though  usually  non-linear  differential  equation  for  the 
position  of  the  interface  Z(t),  that  is  for  the  length  of 
the  column,  as  a  function  of  time.  The  integral  solutions 
vary  in  detail  depending  on  the  conditions  (vertical  or 
horisontal,  flow,  etc.)  and  on  the  assuaiptions  made  with 
respect  to  the  forces  involved.  Newton*  s  law  of  force  applied 
to  the  moving  column  of  fluid  yields  in  general  an  initial 


value  situation: 

p(z,z’  ,Z")  -  0 


2(0)  -  2^ 

Z’  (0)  -  2'  -  V, 


(1) 

(2a) 

(2to) 


(3) 


viher<a  z’  ,  z"  are  the  velocity  and  acceleration.  A  recent 
application  la  that  of  Siegel  QQ  vitio  studied  spontaneous 
capillary  rise  In  reduced  and  aero  gravity  fields  (l.e. 
for  Boute  aiid>lent  conditions  in  space  flight) .  It  should 
be  noted  that  although  rise  %«as  initiated  spontaneously  by 
a  sudden  reduction  in  the  'weight*  of  the  liguidj  the  initial 
conditicns  were  well  defined:  viz.  the  liquid  was  first  per¬ 
mitted  to  rise  to  its  equilibrium  height  under  normal 
gravity  and  cosae  to  rest  {V^  =  0) , 

W.  E.  Brittin  [[731  among  others  analyses  the  dynamics 
preceeding  Siegel*  s  chosen  initial  situation;  that  is,  one 
eiu3  of  a  capillary  filled  with  air  or  vapor  la  touched  to  a 
reservoir  of  wetting  liquid  so  that  flow  into  the  tube  is 
initiated  spontaneously.  Here,  the  flow  initiation  phase  is 
an  implusive  transient  during  which  molecular  forces  of  adhesion 
come  into  play,  ^is  would  be  difficult  to  analyse  especially 
because  of  the  entry  turbulence.  Obviously,  one  should 
attempt  to  set  physically  reasonable  and  mathematically 
consistent  initial  conditiona.  The  present  writers  in  con¬ 
sidering  the  integration  of  variants  of  eqiiation  (1)  by  power 
series,  munerical  methods  and  by  use  of  an  electric  analogue 
device  noted  that  and  z'  ^  as  set  hy  Brittin  led  to  physical 
and  mathematical  inconsistencies.  The  sources  of  these  are 
discussed  in  the  Appendix.  An  alternate  general  p»>cedure 
was  developed,  therefore,  and  the  results  are  presented  below 
for  an  intcgrable  case  of  equation  (1) .  The  duration  of  the 
trauisle-tt  calculated,  as  well  as  the  maximum  velocity 
attained  ivr  all  possible  values  cf  the  viscosity,  density 
and  surface  tension,  dependent  on  the  initial  conditions 


(4) 


subject  only  to  the  linitstlons  of  any  psrtieulsr  nodsl 
(eq.  1) .  In  consequence  an  upper  bound  for  the  naxlnun 
velocity  Is  derived. 

INITIAL  CONDITIONS  AND  THE  TRANSIENT 


A  definite  tine  zero  la  not  obvious  when  one  treats 
the  Initial  liquid  entry  Into  a  (horizontal  or  vertical) 
tube.  Let  US  assune  that  the  contact  between  tube  end  and 
tlie  liquid  is  effected  under  near  equilibrium  conditions^  l.e. 
quaslstatlcally.  Specifically,  let  the  velocity  of  the  tvd>e 
relative  to  the  reservoir  of  liquid  be  zero  the  Instant 
capillary  flow  Is  established.  Then  the  Instant  tdten  notion 
of  a  column  of  liquid  laqpends  would  be  an  appropriate  and 
unaniblguous  tine  zero  reference  point,  and  so  by  definition: 


V 


(3a) 


However,  before  estimating  Z  ,  w  should  note  the  dual  aspect 

o 

of  equation  (1) .  It  Is  a  Xlnematlc  equation  for  the  position 
of  the  advancing  Interface,  but  (more  basically)  It  Is  an 
equation  of  dynamics  for  the  mass  (In  units  of  p)f  r^)*  An 
effective  Inertia  must  be  ovezcome  when  motion  smarts,  of 
magnitude  at  least  that  of  the  mass  of  the  vapor  in  the  capillary 
This  quantity  Is  a  minimum  for  the  Initial  mass. 
ConveAion  to  the  length  dimension  Z,  which  refers  to  the 
liquid  gives  Z^zs  to  denote  that  part  of  contributed 

by  the  vapor.  It  is  a  vbxy  small  quantity  but  aliMtys  larger 


than  zero.  No  Z. 


is  a  weak  lower  bound  for  Z_  since  the 

o 


effective  mass  includes  liquid  above  and  below  the  entry  end 
of  the  capillary.  A  better  minimum  value  is  calculated  by 
including  the  mass  of  liquid  above  the  bottom  of  a  meniscus 
with  volume,  say,  1/3  /f  r^,  and  by  equating  this  to  a  cylinder 


(5) 


of  height  zj*  with  voXum  Iff  z^zj"  tzom  which  -  r/3. 

O  0^0 

ConsoquentXy,  the  totaX  initial  mass  ”  )fv^(PG  *  *■ 

Tfr^p  from  «hichi 

Z  •  +  2  ^  •  Pg  ''^^PL  (3b) 

o  o  o 

F  «  0  and  the  conditions  (3a«b)  oust  predict  an  Impulsive 
transient  of  very  short  duration  in  agreement  with  experience. 
Initially,  there  should  be  a  brief  period  of  positive  acceleration 
tiihlle  the  velocity  attains  Its  maximum  value  at  ^Ich  time 
Z"  a  0.  This  marks  the  end  of  the  Initial  phase.  Thereafter 
retardation  should  predominate  with  Z*^  0,  ultimately  Increasing 
to  zero  at  tlme^lnflnlty. 

To  check  the  dynamics  let  us  choose  a  simple  case,  a 
horizontal  tube  with  the  fluid  column  subject  only  to  the 
capillary  driving  force  and  the  viscous  resistance.  Then  the 
time  rate  change  of  momentum 

d/dt()fr^  p  2Z  )  -  2)frocos  9-  Qffii  ZZ* 
vhere  ^  is  the  surface  tension,  0  the  liquid-solid  contact 
angle.  The  0  capillary  driving  force 

with  9  assumed  constant,  and  QulfZZ  is  the  Poiseuille 
steady-state  viscous  drag.  (The  subscripts  have  been  dropped, 
since  all  quantities  now  refer  to  the  liquid).  Both  force 
terms  are  approximations  since  9  is  known  to  be  velocity 
dependent  [[llj  acceleration  terms  of  the  viscous 

resistance  have  been  neglected  ^l]]  .  Turbulence  at  the  entry 
end  may  be  isqportant  in  the  initial  stages.  These  factors 
would  introduce  additional  retardation  into  equation  (4a) . 

Equation  (4a)  takes  the  form: 

(Zz' ) '  +  Qzz'  -P-0  (4b) 

with  P  -2o  0  ^  convenient  to 
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I 

I 

I 

I 

I 


introduca  tha  dlaanaloolaaa  variablas  T  and  N  defined  by 
t  *  T/Q  and  W  «  ZB/'iT?*  Equation  (4b)  then  becomes 

(ww  )  +  ww'  -  1  a  b  (4c) 


ahich  Integrates  ianediately 

W^/  2  *»  T  •  C.  e*^  +  C_  (Sa) 

^  •  I 

The  Initial  conditions  w  and  W  detaxnlne  the  integration 

o  o 


constants  so  that 

W^/  2  “  T  +  (1-W  w' )  {  e"®  -1)  +  /  2 

o  o  o 


(5b) 


v.^ich  Is  a  more  useful  fox*m  for  discussion. 


The  initial  states  of  the  system  may  now  be  examined 

I  V 

for  a  0  and  for  a  range  of  values  ^  0.  may 

be  as  small  as  we  please*  and  the  behavior  of  the  solution  as 
becomes  very  snail  is  of  particular  Interest  for  our  Initial 
value  problem.  Figure  1  shows  the  dynamics  charactcuristlc 
for  ^  0  ('^<8  particular  case  Is  1)  •  acceleration 

W*  Is  initially  (1/W^)  as  the  lav  of  motion  requires*  axuS  then 
declines  to  zero  marking  the  end  of  the  initial  stage  (T^,) . 

As  noted*  It  Is  negative  tiiereafter.  The  Initial  velocity  Is 
of  course  zero*  then  rises  to  a  maucinum*  at  the  end 

of  the  initial  stage  and  thereafter  decays  going  to  zero  at 
time-infinity.  The  mass  of  the  system  Is  Initially  »  1 
while  the  advance  of  the  Interface  W  »  W«4V^’*0 ;  axMl  both 
steadily  increase  with  points  of  Inflection  at  the  end  of  the 
Initial  stage  vd^ere  KT^). 

Of  interest  Is  the  dependence  a)  of  *  b)  of  the 

RAXe 

duration  T  of  the  build-up  phase  and  c)  of  the  jnerement  of 

ID 

mass  (or  advance  of  the  Interface)  d  W  »  W(T  )-K^  on  our 

m  no 

parameter  w  .  First*  the  diuratlon  T  Is  given  as  function  of 

O  ID 

by  the  equation  arising  from  the  condition  W  »  0: 


(7) 


-  <  2  *, 


m 


+  -1«0 


<6) 


Then  V  (w  )  and  ^  If  (W  )  can  ba  calculated  froa  equation  (Sb)  • 
no  n  ^ 

The  three  functions  are  plotted  In  Plqure  2.  There  Is  always 
a  bulld-up  stage  since  H*(T)  >■  0  has  a  positive  solution  which 


0  as  H  — ^ 
o,  ' 


0.  In  the  llnlt  the  actual  velocity  is 


1/2 


|(^  cose)/(pz^^' %  a  definite  upper  bound  defined  by  equation 

(4a).  This  bound  Is  the  sasM  for  vertical  rise. 

The  drag  tern  due  to  the  turbulsnss  has  been  given  as 

%(Z  )^  and  as  proportional  to#  **^(Z  ^  analysis 

of  Langh&ar  ClO^  of  a  sinllar  problem  leads  to  .64  (  Z"  )^. 

In  general  for  any  situation  one  can  write  (In  the  W,  T 

»  2 

variables)  G(T)  (W  )  with  0(T)  some  positive  fvmctlon  and 
0  (T)  »  0.  This  tern  destroys  the  linearity  of  equation  (4a) 
but  must  lead  to  a  snaller  value  for  V(7(|f)  "  V  max.  One 
result  however «  follows  Isnediately!  it  can  be  shown  that 
the  upper  bound  for  V  max.  l.e.  the  limiting  value  as  0 

now  becomes: 


pr 


1/2 

T-0 


that  is,  it  becomes  smaller  by  the  factor  Q/  ^1  <«(0)« 

A  representative  system  Is  a  capillary  of  10  ^  raAus 

.  3 

contalnixig  a  liquid  of  0.12S  poise  viscosity,  1  gram/cn 

For  these  parameters 


density  ^uld  50  dynes/cn.  surface  tension, 
z  W/10*,  t  »  T/10^.  For  Z  -  r  -  10”*cm.  H  -  1. 


The  bulld-xip  tine  T 


-  r  -  10”*cm.  ^ 

o  o 

0.14  milliseconds  and  the  maximum 


velocity  is  50  cm. /sec.  The  advance  A  W  during  the  period 

ID 

Is  0.8W  ,  l.e.  0.8r  cm.  These  figures  have  reasonable 

o 

orders  of  magnitude.  One  should  bear  la  mind  the  limitations 
of  the  model  (Equation  4),  and  the  ideal  physical  ccHidltlons 
prescribed.  The  effect  of  whimsical  factors  such  as  the 


(8) 


deviation  from  th«  naar-aquilibriun  conditions  at  tha  flow 
initiation  mtagm,  rugosity  of  tha  tuba  and«  ate.  cannot  ba 
estlaated.  Uhlass  tha  intarast  of  tha  axparinantar  is  in  tha 
initial  stagas,  it  will  ba  convaniant  always  to  bava  sona 
liquid  in  tha  tuba  at  tiaa-saro  as  in  tha  axpar inants  of  Siagal 
1^^  .  fiqposing  this  condition  rsircuavants  tha  analytic 
difficulty  of  treating  tha  initial  davalOFOMunit  of  tha  inter¬ 
face  curvature  and  the  attendant  transient  end-affects  %diich 
occur  when  the  tuba  end  is  first  brought  in  contact  with  a 
liquid  reservoir. 
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APPENDIX 

Brlttin*  8  equation  is 

(zz')'  +b(zz’)  (Z)^  ♦<Z^ciKJ  (7) 

with  -d  the  same  physical  quantity  as  P  in  equation  (4a) . 
the  initial  conditions  are  set  in  this  manner:  ''Z(0)<>0''4 
ii  "Ihe  initial  velocity  is  assumed  to  be  finite  as  it  must 

be  physically!  and  thus  is  obtained  from  the  differcaitial 

k 

ecniat4.on  bv  setting  Z«  0.  (-4/Sd)  ",  (Ihe  italics  and  the 

numaration  and  i^i  by  the  present  writers) . 

We  note:  that  (1)  the  '  £eed~back'  procedure  in  ij^  is 
not  valid!  since  both  conditions  must  be  set  independently; 
that  (2)  Z^  a  0  sets  the  initial  mass  and  the  initial  kinetic 
eneargy  equal  to  zero;  and  that  (3)  Z  »  0  suggests  the  initiation 

I  1*  ® 

of  flow^  vdille  =*(~4/5  d)  ^  which  is  about  80  cm.  per  second 
in  the  representative  system  mentioned  above>  really  applies 
to  a  later  stage  of  the  notion. 

To  illustrate  this  difficulty  we  proceed  as  foilotra: 
equation  (5b) ,  which  is  the  general  solution  of  (4b)!  Includes 
all  pairs  of  initial  conditions  and  not  only  those  applicable 

t 

to  spontaneous  Inibibition.  Let  vts  wlte  W  »  ^  ^ 

o  o  , 

thus  define  a  function  9  of  three  variables.  Kow  take 
arbitrary  but  bounded^  and  exanlzM  the  behavior  of  ^9/^T  near 

t 

(0,Q,vi^  ).  First  let  T  go  to  sero!  yielding: 

Lime  *W^* 

W^->  0  ^  T  T»0  ® 

as  we  wish  i.e.(^9/"^T  at  (0!0,l#  )  ■  .  Also  W  »  l/H 

00  00 

:4)ich  id  JKewton’ s  law  of  motion  in  dimensionless  units  for 
unit  force.  Now  let  0.  IhJ.s  yields  (there  is  a  rcncvabla 

indeterminacy) : 

Lim 

0 


Thus  the  limit  depends 
Is  excluded,  the  limit 
path,  and  consequently 
behavioE  at  this  point 
togral  of  equation  (7) 
expansion  about  t  =»  0, 


on  the  path.  However,  if  W  *0 
of  1  is  w^lndependent  of  the 

(0,0, )  Is  a  singularity.  The  mis 
is  immediately  apparent  for  the  in* 

I 

near  (0,0, ):  the  power  series 
contains  terms  with  1/Z  as  a  factor 
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FIGURES 

Figure  1.  Initial  stages  of  capillary  rise.  Plot  of  Inter¬ 
face  advance velocity  w  and  acceleration 
w* '  ys.  T.  Initial  xoaas,  1  unit. 

Figure  2.  Dyii3sjntcs  of  initial  laipulnive  transient  as  a  function 

of  IT  j  the  effective  initial  mans  (or  length  of 
o  , 

coluion) .  Plot  of  naxlmuin  velocity  attained  W 
(max)  >  T^^tiroe  to  attain  W  (max)  ^  and^W^  the 
interface  advance  vs  W  (log  scale) . 


Figure  1 


10-3  10r2  10-A 


Sp«6tr«X  Aaalysifl  of  Stxoeturo  of  fbxouo  liidlo 

By 

Borxy  Fora 


aiwflwtawii 

Hio  gooMtrloal  goantltioa  gaaonlly  aaod  to  oharaetarlao 
a  poxoua  aMdiiaij  aoeh  aa  poroaity  and  ap*eifle  muzfmam,  aro 
ovar«all  avaragaa  and  ara  not  auffLeiant  to  daacribo  it 
exactly  or  identify  it  uniquely.  Zt  would  be  de8iand>le  to 
have  criteria  by  which  it  can  be  decided  by  a  aiaqple  ahalyBls 
of  photoaiicrographa  whether  two  porous  nedia  are  alike  or  aot» 
and  to  what  degree,  teaie  prellainary  ideas  in  this  direetion 
ifere  outlined  by  Fare  and  d^ieidegger^^^ . 

Zt  is  obvious  that  in  order  to  define  *lower<-order* 
structure  factors  (related  (say)  to  poxooi^)«  recourse  oust 
still  be  BMde  to  statistics.  Ns  first  coostruet  a  function 
which  records  not  only  the  porosity  but  the  "lower-order* 
density  factors  as  well.  Aa  arbitrary  line  is  wssuned  to  be 
draon  through  a  given  porous  nediuai  whose  gsosietry  is  to  be 
described.  Points  pn  the  line  are  to  be  defined  by>giviag 
their  are  lengtit  1  froo  an  arbitrary  Chosen  origin.  ..  ihso^ 
for  eertain  values  of  •  the  line  villv  pose  thicongl^.. void 
spaces;  for  other  values  of  •  it  will  pass  through  filled 
spaces.  Ns  tbrn  introduce  a  fuaetion  f(B)  ef  the  aare  length 
$  defined  as  followet  the  value  of  f  is  defined  as  tl  if  the 
line  at  the  point  passes  through  a  void  ppaoej  it  is  defined 
as  equal  to  -1  if  the  line  at  8  passes  through  filled  iqpaoe. 
tte  fuaetion  is  a  raadon  "rectangular  wove*.  (Bee  f  in  figure  1)  • 


w 


Wo  note  the  tisean  t{9}  i.»  related  to  the  porosity  P  by: 

f  2P-1  .  (1) 

ioer  in  the  Xiasi  t  the  intearval  over  ^ich  f  Is  defined 
is  increased;  the  ^^lAtionship  (1)  vlll  texid  toward  equality. 
<!h&  definition  of  ^  can  of  course  b/'.  jocx^ified  so  that  the 
jEtSKi  'woulf  yielt'  t*  tijcactxy  biit  there  i  v  no  advantage  for 
lifter  analysis. 


the  first  te?a!«  of  a  Fourier  eispanaion  of  any  function 
yields  its  Bteeui  value,  &o  that  (1)  svtggests  a  spectral 

i/e  write  fozsally 


dccoaqpositlon^^  f  (S^ . 
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'imere  1  is  the  interval  analysed  and  h  B  are  the  Fcorier 

Xh  K& 

co»efflcients  to  d^^ter^ned*  Uiey  are  gi  risn  by: 

pL  <  ^ 

i  \  ^  2  m  dB  j  s  ®  il  i 

^4  ^ 


£is)jQAM,tlJBSL 
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Use  sets  of  co-^efsicienfes  (A  ,fi  )  n*K>,  1,2,3, ....  these 

II  n 

characterise  the  porouis  ):«^iuai.  It  la  clear  that  any 
differences  in  suifstruc  tore  between  t'*eo  pon:rus  media  will 

be  reflected  in  dlff ericas  in  some  of  the  sets  CA  ,B  ). 

■  — '  n  o 

For  exaaple,  if  the  porosities  P  are  identical  then  fox  n°C, 
we  will  have  A^»2?.>1,  boti;  cases, .  but  other 

sets  will  differ  in  ^>jTr<^.ral.  ^e  emslysis  can  be  siapllfitid 
conoepUially  by  the  ,r;«.^llowixig  conaiderationds  the  terms 
in  (3)  can  be  coBbir/.e»f  in  pairs  into  terms  containing  a 
phase  anglej^uhich^ho^v^'  ex  df^pend#  cn^  choice  of  orl^'in, 
and.  thus  is  not  ro  our  anaXv'^.l^o  new  co-efficients 


are  given  fayt 


(3) 


(3) 


TlM  C  th«a,  «r<t  th«  tru«  mpt^trml  A«n«itlM*  (Zt  ftolZowi 

la 

tiMit  c  ■  2f*l.  )  . 
o 


l!b«  ftiactloa  f  (s)  iM  only  ooattoMoiui.  Mom 

it  would  bo  ooofttl  if  tbo  C  r«f  l«ct«d  tho  oubotrustuxo 

A 

ia  M  slB^lo  a  mdcy  am  posaiblo,  vis.,  hy  th«  seat  rapid 
conrargbuce  of  tbo  sarioa  {2),  ao  Ibat  tho  C  fall  off  aa 

a 

rapidly  a«  poaslbla  with  iskcraaal^  n.  Iba  rata  of  oonvaxgaaoa 
ia  blQbar  for  contlnuotta  fuBctloaa.  Botindiaf  off  tba  ooraara 
of  f  (8)  in  soaa  cowaiatant  faahioa  aaggaata  itaalf.  aacli 
aquara  wava  ia  f(0)  waa  raplaead  by  a  aiaa  wava  of  tte  aaaa 
araa  and  of  tba  aana  baXf«wava  langth.  fba  aaw  funetloa 
p^i)  ia  alao  abowa  in  figure  (I) .  iba  diacoatiaatiaa;  ara  aov 
ia  tba  first  derivative.  A  siapla  aaalyaia  ahowa  however* 
that  the  sine  wave  arcs  introduce  axtranaeua  **aoiae". 

fha  f  and  F  of  figure  (1)  are  tba  atrueture  fuactiona  of 
a  aandatona  of  26.3%  porosi^.  ibe  line  ia  vwy  abort  (0.17eai). 
The  firat  14  valuaa  of  c  for  F  and  F  wera  oowpetad  aod  are 

A 

abovn  in  table  (1)  •  «be  rectangular  wave  daaaitiea  are  plotted 
in  figure  2. 

Figure  3  auggeata  bow  tba  reaulta  of  tlM  apeetral  aaalyaia 
can  be  preaantad  gri^ically  to  provide  a  atrueture  "pcofila*  of 
a  poiDua  sMdiue.  Of  particular  interest  would  be  tba  peeks 
in  the  c.  if  amy,  Zt  abould  be  noted  that  tba  langth  uaed 

A 


(4) 

ill  this  particular  eaaa»  (O.ITcb)  iaivolulaf  a^p  12  altamationa 
ot  void  and  fillad  ragioaa  la  proMblp  tao  aaall  to  mUeo 
fiigura  a  aeeurata,  altioagb  tba  caXculatad  poxoaitp  ia 
28.8%  in  oeovartaoB  wilii  tta  kaeua  ualua  of  88.3%.  Oaa  Mp 
oamarfo  on  tba  pcafila  by  aaalyaiap  attoooaalra  laagar  lauptlia 
Hm  ipaetral  daaaitiaa  mut  atom  bo  plattad  uo.  tba 
uaua  auii>ay  V  **  'tiMt  tba  uariona  paofilaa  can  ba 

ahaim  on  ana  grtpb.  tba  dagraa  rata  of  convargianea  can  ba 
noted  at  aacb  ataga. 


Slaea  tiia  funetlaa  f  (•)  la  alttar  ‘*'1  or  •>!»  eaavutatleii 


of  ttia 


and  la  rolativaly  alapla. 

*  of  tW  witb  tba  dlaoentlimltlaa  at  a^  b  . 

o  o 


a,  b_f»**a. 


•1  '1* 


■n-1 


for  tba  atrlpa  Ca  b_) » 

o  o 


(a^b^)  •  •  * 


(a^b^)  aat  f(«)  «  +1;  for  tba  atrlpa  Cb^*j^)  (b^a^)..* 
aat  f(8)  «  --1. 


(11  taho  eoLa  alapla  fdi«t- 

“1  -#*^1 


(5) 


fflth  a  alnllar  fosanla  far  tba  a^.  foraala  (•)  eaa  ba  fortiMr 

m 

radaead  .for  aiuaarlcal  warka  alaoa  tba  Intavgratioaa  can  ba  par- 


tba  oaXcalatlana  can  ba  progri— ad  galekly  far  a  dlgltlal  aaapatar 
far  tba  ganaral  cmoo;  la  aacb  tnatanca#  a  Hat  of  idm 
ublab  aaa  ba  road  off  a  pbatnal  nmgraib,  aad  tba  laagtb  1^ 
ana  Id  ba  tba  only  lapata  aaeaaaary. 


(1)  faro  and  Icbaldaggar  *8tatlatleal  Oaonatxy  of  foxata  Madia* > 
J.  oaa.  baa  S§,  (10)  1981. 
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mu  1.  Usmotxml  Oo*^floi«nta  of  a 
Foreoa  NadluB 


Fuaetlon  £ 

rimetlon  W 

a 

C 

c 

a 

0 

0.464 

0.465 

1 

0.267 

0.268 

2 

0.606 

0.714 

3 

0.272 

0.457 

4 

':*.i54 

0.052 

5 

0.277 

0.289 

6 

0.232 

0.498 

7 

0.244 

0.474 

6 

0.439 

0.478 

9 

0.260 

0.247 

xo 

0.3S3 

0.410 

11 

0.142 

0.110 

12 

0.158 

0.154 

13 

0.229 

0.278 

14 

0.088 

0.085 
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■OMBRZCAL  ZMTBOIATXOIII  OT  8T0Klt-mVZIR  IQOIkTZOn 

»ar 

Harry  Fara,  O.  Oalgar  and  w.  Haaa 


nia  poaitlon  la  oftaa  takaa*  that  tha  aqaatlooa  of 
stokaa-Haviar  ara  parfaetly  gaaaral  for  tha  daaeriptlon  of 
all  typaa  of  flow  (a.g.  tuxbolant  aa  wall  aa  laaiaar)  la 
May  apaea  (a.g.  porooa  aadia  aa  wall  aa  unbouadad  apaca)« 
and  vmdar  all  eonditioaa (nnataady  aa  wall  aa  ataady  atataa) . 
lha  ganarality  of  thla  propoaltioa  (of.  S^lichtlng,  loundarv 
Lavar  lhaorv.  Iiearaw-Bill»  IMO)  ia  Ita  own  llaitation, 
iaaanooh  aa  iatagrala  of  tha  raaultant  aonliaaar,  high-ozdar 
partial  diffaraatlal  aqnationa  of  notion  ara  not  known  In 
cloaad  fom  avan  If  wall  oondltlonad  Initial  and  houadary 
condition  could  ba  apaelf lad  (which#  for  tha  gaaaral  caaa# 
la  aot  ppaalbla).  On  tha  othar  hand#  approxlnatlona  aoaiatlnaa 
ara  poaalbla  which  load  to  a  fair  daacrlptloa  of  tha  dyaanica 
of  haalcally  eonpian  fluid  aotlono  (a.g.  lha  Oaaaa  approxlnatlon# 
tha  Praadtl  "bouadary-layar”  prooadura#  ftokaa  *oraap  flow^# 
ate.).  Wa  ebooaa  now  to  oonaldor  aaothar  poaalblll^#  aaawly 
tha  "bruta-forca”  aunarleal  latagratlon  of  tha  itokaaHlawlar 
aqaatlooa  by  raaort  to  flalto-dlffaraaea  aaaloga  aa  awaluatad 
by  tha  uaa  of  hlgh«apaad  oonputar  proeadoraa. 

Aa  aa  alanontary#  alhalt  haorlatlcally  atlanlatlag  eaaa 
to  ooaaldar#  wo  chooaa  tta  pcoblaai  of  Figaro  1#  wharo  a  aat 
of  oooplad  platona  (aaparatad  by  a  dlataaea#  d  nova  ia  ataady 
notion  through  a  oyiladrleal  tuba  of  radlaa#  a.  Lot  aa  aiqr 
that  tha  oontalaad  fluid  la  inooapraaalbla#  and  that  tha 
valoclty  of  platen  nowiant  la  fhia  eaaa#  which  ia  naha* 
lagfal  for  aa  applicatloa  to  tha  foantaia-affaet  prohlan 

I  Ilhll.M  III  ■■  I  1.1. 1  ■Kll 

*of.' Blrkheff#  hydrodyaaBiea# 
faooad  ■ditioo# 


(2) 


po««d  by  ]lo««  (Vatur*»  Jbly  15»  I9€l),  Is  rspMssntsd  byt 


•ad 


^  ti9/Dt  •  “  'V  ^ 

dlv  ^  -  0 


(1) 


•bleh  •qiiatlons  rafar  of  eourao  to  •  fixod  Balorlaa  eoordiaato 
•ystMi.  Xa  Bquation  1«  v  is  tho  voloeity  vootor  hsviag  tbo 
axial  coaqpoaaat,  u,  aad  tba  radial  coaqponaat*  w;  aad^ 

aro  thm  fluid  dansity  and  viscosity,  raspaetivaly;  aad  P  is 
thm  i^assura. 

Taking  tbasa  aquations  in  cylindrical  eoordinatas, 

(r,  X,  ^  ),  and  fixed  to  tba  eouplad  aoving  pistons;  aad 
furtbar,  assuaing^  indapandant  of  ^  and  tiaa,  aa  obtaiat 

•r®x  “  ®x*R  ■**  “x  *  *  P 


is  tba  oparator  giving  D^8  «  8  ■  T 

(i.a.  O^T  is  tba  bibaraonic  of  sf« 

8*-  ru 

8  *  •r  w 

: .  .V. 

aad  aliara  tike  subscripts,  R,  r  and  x,  raprasaat  diffaraatiation 
of  the  iadleatad  function  altb  raspaet  to  tba  subaoript. 

Tba  no-slip  boundary  conditions  for  tba  casa  of  viscous 
flow  (cf.  Figura  1)  ara: 

8  *0  for  R  *  a®/4  aad  for  x  ■  dtd/f 

X 

8^  •  0  for  X  -  1  d/2 
8_  •  for  R  -  0^/4 

With  tbasa  boundary  conditions  iaposad,  tba  staady-stata 
solution,  8  •  8(x,r},  is  to  ba  found  by  a  auaarioal  itarativa 
proeadura  with  tbm  finita-diffaraaea  analog  of  Rquatioa  2. 


(3) 


A  prttXlBliiary  study  iadiostss  that  tbs  Itsrattons  aaeassary 
to  ol»tain  oonvaxganea  can  Im  uadartakan  with  aeonosqr  oa 
aadlUB  spaad  ooaputors.  Plgura  2  daplcts  sehsaatioally  tha 
flow  pattara  whioli  Is  stpgMtad  hf  physical  iatuitioa.  To 
ba  BOtad  Is  tha  fact  that  sysaatrias  otibar  than  axial  syapstcy 
cannot  ba  pradlctad  ariorii  aoraevar*  it  is  likaly  that  ttia 
sinpularitias  at  tha  piston  pariaatars  (arisiag  baoausa  idia 
valocity  coaponantf  u,  is  not  siagla<^aluad  at  thasa  points) 
will  introduca  a  difficulty. 


(4) 


HBTMORK  JUfliLOOOBS  FOR  MOdUBIB  OF  FLOW 
TBSOOaa  FOR0D8  MBOZA 

IZZ  note  on  a  "DiocrlBlnating"  Routine  For  A  Step-Wlaa 
Daaaturatlon  Frocaaa 

By 

Barry  Fara 

Ilia  applicability  of  alactric  nataork  analoguaa  to  the 

study  of  flow  through  porous  aadia  is  lisiited  by  the  large 

nusdbar  of  cos^nants  required  to  sensibly  represent  an  actual 

sasiple  of  porous  siaterial.  This  is  true  of  course  for  physical 

(1) 

models  constructed  with  resisters,  etc.  but  applies 
as  well  to  "networks”  %diicb  are  to  be  solved  by  pencil-end- 
paper  <q>erations.  we  are  ooocemed  with  the  latter  case, 
for  which  this  limitation  can  be  overcome  in  part  by  the  use 
of  high-speed  digital  computers.  Ihe  flow  chart  of  a  programme 
for  solving  a  steady-state  single  phase  problem,  written  for 
the  universi^  of  Illinois  computer  "Zlliac”  was  described  in 
section  Z.  Zn  section  ZZ  (Seventh  Ouarterly  Progress  Report) 
an  analysis  was  presented  of  a  procedure  suitable  for  a 
stepwise  desaturation  process  i.e.  when  a  wetting  fluid  B 
replaces  a  liquid  A. 

A  critical  operation  in  the  procedure,  whether 
effected  manually  or  with  a  computer,  is  keeping  track  of 
possible  paths  for  fluid  motion.  At  the  i^  stagei^(say}  the 
"pressure”  isqposed  is  such  that  (1)  pores  of  sices  *■  S^  are 
now  possible  flow  channels  for  tiie  invading  fluid  (see  ZZ) ; 
b>it  for  B  to  replace  A  in  a  given  pore  we  must  also  have  (2) 
an  exit  path  to  the  sink  boundary  for  A  and  (i)  a  continuous 
path  of  pores  of  sices  ^  S^  from  the  source  boundary  to  the 
given  pore.  Zf  these  three  conditions  hold,  then  a  pore  is 


1. 

2. 

3. 


allowad  to  "pop*  !.*,  fluid  P  io  coiuildoffod  to  havo  ro* 
placod  fluid  A.  Application  of  thaao  criteria  is  sisqpla 
but  tedious  in  nanual  eork  *  •  For  Mchine  processing  a 

sequenee  of  disclaiaating  binary  operations  resulting  in 
yes-«o  decisions  is  required  ^ieh  will  duplicate  tae  results 
of  va.sual  discrininations. 

Consider  Figure  1,  a  Miniature  rectangular  network  con¬ 
sisting  of  resistors  (conductances)  represented  by  dashed 
and  solid  lines  and  linked  at  the  nodesj  (the  circles) .  The 
solid  lines  represent  pores  now  containing  fluid  the  dis¬ 
placing  phase*  idiile  the  dashed  lines  r^resent  pores  still 
containing  fluid  A.  The  source  boundairy  is  at  x>0.  The 
conductances  to  the  lift  of  (solid  lines)  %dtich  are  not 
part  of  the  network  conveniently  fix  the  line  as  the  source 
of  fluid  B.  The  desaturation  process  is  at  soms 
intemediate  stage  (say  the  i^)  with  the  network  partitioned 
as  shown.  The  conductivities  of  each  can  now  be  calculated 
(see  X) . 

We  now  assune  that  the  pressure  is  increased  so  that 
pores  of  sise  ^  8^  are  now  possible  channels  fOr  B;  these  stay 
include  sone  of  sises  ^^^^  which  did  not  *pop”  at  earlier 
stages  l>ecause  conditions  (2)  and/or  (3)  were  not  satisfied. 

The  scanning  procedure  begins  at  this  point.  It  is  clear 
that  we  siay  not  proceed  in  any  arbitrary  way.  It  is  the 
logical  (snehine-prograasaad)  attack  of  ^is  problen  which  is 
to  fom  our  future  work. 

&il£SCSB£ttL 
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Contact  Anglo  Aeooloratlon  Oopondoney 

By 

N.  Chaudhari  and  waitor  Soao 

Ao  a  eoatinuation  of  tho  work  doooarlbod  hy  tho  oritar 
la  tho  sixth  and  oovonth  Qoartorly  Pxogrooo  Aoporto  of  this 
aorlos,  ojq^lMnto  ooro  undortakoa  ohoro  tho  dooay  of  m 
aanonotor  wao  oboorvod*  Zn  ouch  a  aymtm,  aceoloratloao- 
docoloratlono  occur^  so  that  tho  ohsorvatlon  of  tho  Intor- 
faco  position  of  tho  aanonotrlc  fluid  as  a  function  of  tlsM 
provides  tho  data  against  i«hlCh  hypothosls  can  ho  tostod 
about  tho  dopondency  of  tho  contact  anglo  of  aovlng  fluld^fluld 
Intorfacos  on  tho  voloclty*  and  on  tho  tins ■ dorlvatlvo  of  tho 
voloclty  of  Intorfaeo  sovoaent. 

Zn  tho  oaqporlsMnts  roportod  horo«  data  ooro  Obtalaod 
using  aanoswtors  ohoro  tho  dlaswtor  of  tho  tdbo  anas  olthor 
%fas  0.185  cm,,  0.200  ca.,  or  0.200  ca.  Dow  Coming  (VaalMr  200) 
silicon  oils  aoro  usod  as  tho  aanoastrlc  fluid*  having 
vlscosltlos  olthor  of  10*  7*  5  or  3  omtipolsos.  iha  abolo 
systoa  aas  kiqpt  undor  tsaporatur^  control  at  25*  C.  Tho  position 
of  tho  trailing  aanoastor  Intorfaeo*  Z*  as  a  function  of  tlao* 
aas  rocordod  with  a  Paths  (liiho)  16  sai.  aovlo  caaora. 

zn  all  casos*  tho  oaqporlaontal  data  could  ho  r^prosontod 
by  a  llnoar  rolatlonOhlp  hotaoon  tho  logarltfaa  of  tho  ratio 
(I^Z^-Z)  and  tiao*  ahoro  aarks  tho  tlao  infinity  position 
of  asnnaotsr  balaaeo.  According  to  theory*  dovlatlons  frea 
llnoarlty  only  aro  to  ho  saqpootod  during  tho  first  allll* 
sooonds  of  tho  iapulslvo  transient  ahlch  folloas  tho  initiation 
of  notion,  ihoory  also  prodlcts  that  tho  slops  of  tho  straight 


(2) 


line  function  will  liavo  •  high  value  if  the  advancing 
and  receding  contact  angles  in  the  leading  and  trailing 
nanoaeter  ams  have  the  saMi  value  through  tte  decay;  vbile 
a  sehifibly  lower  value  for  tiie  slppe  is  predicted  on  the 
assumption  that  the  receding  contact  angle  has  sons  fixed 
(tiae^invarifuat)  value  (say,  zero),  and  on  the  assuaqption 
that  the  advancing  contact  angle  has  a  tine^dependency  repre* 
sented  by: 

cos  ^  ^  •  1  +  +  MjZ*  (1) 

where  ^  ^  is  the  advancing  contact  angle  of  the  leading  in* 
terface,  and  z*  and  Z”  are  tiie  velocity  and  the  acceleration 
(deceleration)  of  the  notion,  and  and  are  constants. 

iguation  1  is  arrived  at  by  postulating  that  ^ 
both  veloci^  and  acceleration-dependent,  as  given  by  talcing 
the  first  two  terns  of  a  double-Taylor  series  expansion. 

This  fbm  was  suggested  by  taking  note  of  the  observatlmis 
of  Bose  and  Heins  (cf .  previous  progress  r^^rts  of  this 
series,  cf.  also,  J.  Colloid  Science,  Pecenher,  1961)  that 
the  cosine  of  is  linearly  related  to  velocity  for  steady- 

state  notions. 

Xn  any  case,  a  neasurenent  of  the  slope  of  the  eiqperi- 
nental  data  plot  of  In(Zj^Z^-Z)  versus  tine,  according  to 
theoary  (i.e.  a  closed  analytic  solution  of  the  nanoneter 
problen  is  known)*,  provides  a  way  to  evaluate  the  constants 
of  Bquatiem  1.  Table  Z  shows  results  obtained  for  the  3 
esntistoke  oil,  where  the  values  of  and  l^ve  been  obtained 
by  a  least-squares  fit  of  the  data  nakiag  use  of  an  ZBN  650 
progran.  Here  it  is  seen  that  the  nagnitude  of  the  acceleration- 
dependency  (Mj)  is  of  the  order  of  only  lO-peroent  of  the 


*cf .  Ml  thesis  of  H.  Chaudhari,  oniversi^  of  Illinois, 
Psbmary,  1962. 


(3) 


Velocity  dapMadwiey  (N^)«  and  ia  of  oppoalta  sign.  Vha 
HMaaing  of  thla  ia  that  tho  valocity«dapandancy  (aa  naaaurad 
by  nagativa  valuaa  of  M^)  aarvaa  to  alow  down  tha  nanonatar 
daeay  in  acooxdanca  with  tha  principla  of  LaChataliar  aa 
pradictad  by  Soaa  amd  Baina  floe.eit.) ;  whila  tha  aeealaration- 
dapandancy  acta  aiailarly  baeauaa  in  tha  praaant  inatanea  of 
nanonatar  dacay«  tha  accalaration  ia  nagativa  (i.a.  tha  dacal* 
aration  ia  poaitiva)  aftar  tha  firat  iaqpulaiva  inatanta  %dian 
the  notion  ia  initiatad. 

Tabla  X  alao  auggaata  that  both  valocity  and  accalaration- 
dqpandancy  incraaae  aa  tha  tuba  aiaa  incraaaaay.  rapraaanta 

a  trend  which  ia  pradictad  by  tha  argunant  that  tha  capillary 
forcaa  are  laaa  affactiva  in  large  tubaa  (than  anall  thbas) 
to  naintain  high  curvature  of  tha  moving  interface,  tharafora, 
flattening  by  tha  hydrodynanie  forcaa  (i.a.  tha  "fountain  affect* 
of  Rora,  Nature,  July  15,  1^61)  occura. 

Since  Tabla  I  indicataa  that  tha  accalaration-dapandancy 
ia  anall  coa^parad  to  tha  valocity-dapandancy,  tha  neglect 
of  tha  higher  order  and  degree  tama  of  tha  Taylor  aaqpanaion 
of  Equation  1  H'P<ursnbly  ia  Juatifiad,  alnilarly,  there  appaara 
to  be  no  raaaon  to  conaidar  axpanaiona  which  taka  ooa  ^ 
aa  a  function  of  (Z*  ••••)  •  Zndaad,  Ziagal  (J.  Applied 

Nachanica,  June  1961)  haa  even  dropped  tha  tarn  of  Equation 
1  in  hia  analyaia  of  capillary  riaa  (i.a.  a  decelerating  ayatan), 
tharafora,  wa  ahow  in  Tabla  ZX  our  calculationa  of 
(i.a.  tha  valocitv4^andaney)  which  rapraaant  tha  beat 
(laaat-aquaraa)  fit  of  our  data  on  tha  aaawption  that  can 
be  aat  equal  to  aaro.  Bara  again,  a  general  incraaae  ia 
ia  obaarvad  correlating  with  aa  incraaae  ia  tuba  aiaa;  aoraovar, 
alao  appaara  to  incraaae  with  incraaaing  viaooaity.  Both 


(4) 


o£  tt%»m  trends  are  oonsistant  with  thsorsticsl  «xpsotstions 
as  lapliad  above  (c£.  NS  thesis  of  N.  Chaudharl*  loc.cit.) . 


(5) 


Thm  y 


Tdb«  81s« 

O.llC  R. 

—  —  .  B 

O.IC  R. 

B  ..  . 

0.0915  C_  R. 

0.156 

0.0926 

0.084 

0.0122 

0.0094 

0.00036 

TABLE  II 

Tub«  Sixe 

O.llC  R. 

O.IC  R. 

0.0915  C^  R. 

.....  .  _  ....  B  .. 

.  B 

B 

Viscosity  of 

j 

i 

1 

oil 

10 

0.75 

0.94 

SB 

1 

7 

0.5 

0.377 

0.201 

5 

0.65 

0.63 

0.438 

3 

0.28 

0.083 

0.053 

